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ABSTRACT 



Sound production and modification in an Argon discharge plasma 
was investigated theoretically and experimentally. Sound pro- 
duction by ions in the cathode region was shown to be feasible and 
consistent with earlier experimental results. A numerical analysis 
of the theoretical dispersion relation of Ingard and Schulz showed 
that plasma amplification effects will not compensate for acoustic 
losses. The theory of the mobility-limited thermionic diode 
operated as a microphone was expanded. Such microphones were used 
in an experimental investigation of sound production and modifi- 
cation by a discharge plasma. Possible plasma effects on sound 
propagation were observed. The characteristics of sound pro- 
duction by the plasma were found to be in reasonable agreement 
with theory. A departure from theory at low frequencies was inter- 
preted as arising from the effects of particle diffusion, an effect 
not previously observed in acoustics. 
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INTRODUCTION 



I . 

Plasma acoustics deals with the interactions between charged 
particles and sound fields. The problems encountered are different 
from those of ordinary acoustics since charged particles experience 
bulk forces from external fields whereas neutral particles do not. 
When neutral particles are present, they are indirectly influenced 
by the external fields through their interactions with the charged 
particles. The complete mechanism of interaction is very complex, 
involving ionization and recombination as well as elastic and non- 
ionizing inelastic collisions. 

This paper focuses on the problem by studying the effect of 
the interactions on the neutral particles. Two aspects are 
treated. One is the creation of sound by the current modulation 
of a discharge. The other is the modification of the sound^ when 
it propagates through a discharge. The two aspects are closely 
r elated. 

The creation of sound in a gas requires a localized time-de- 
pendent pressure variation. In a weakly ionized gas, the partial 
pressures of the ion and electron fluid components are negligibly 
small with respect to the partial pressure of the neutral gas. For 
this reason, the acoustic wave made in a discharge plasma is con- 
cerned only with the neutral component except when the effect of 



’’Sound” and ’’acoustic signal” or ’’acoustic wave” are used 
synonomously and are tiefined to mean any fluctuation in the 
pressure of the neutr£vl particles. 
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ions and electrons on the neutral component is of interest. The 
pressure, expressed as p = kNT, can be changed by either a change 
in temperature, T, or particle density, N. 

In the positive column of a glow discharge, energy is con- 
tinually fed into the plasma from the electric field, which draws 
its energy from an external power supply. The energy is received 
by the electrons and ions in the plasma as kinetic energy. Since, 
in the aggregate, these particles are traveling at their terminal 
velocity, referred to as drift velocity, the increased kinetic 
energy of the individual electrons and ions is dissapated through 
collisions with the neutral particles. The average rate of energy 
drawn from the field, per unit volume, can be written 

^ 2 

H = H + H. = eN u. + eN.M-.E 
e 1 e e 1 X 

where M* and |x . are the electron and ion mobilities. Since the 
e 1 

electron mobility is much larger than the ion mobility and the 
other terms are equal, the ion term *is negligible. Thus only the 
electrons need to be considered in the energy transfer process. 

In the equilibrium situation, the energy received by the 
electrons must be balanced by losses. The energy distribution of 
the electrons closely approximates a Maxwell Boltzman distribution 
The difference is due to the effect of the electric field and the 
energy-dependent cross sections of the neutral particles with 
which the electrons interact. The electrons transfer energy to 
the neutrals through elastic and inelastic collisions. The energy 
transferred in inelastic collisions is usually transferred out of 
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the system of neutral particles through radiation. The energy 
received by the neutral particles in elastic collisions is thermal 
energy which, at equilibrium is balanced by thermal losses to the 
container and neutral gas outside of the ionized region. This 
energy transfer between the hot electrons and cold neutrals is 
comparable to mixing two gases at different temperatures except 
that the hot gaseous component stays hot. The elastic collision 
losses .ar e. normally much larger than the inelastic losses which 
can be neglected for tliis analysis. 

Thus the processes discussed suggest a mixture of two fluids 
of different temperature which retain their identity in the mixture. 
The hot fluid receives energy from an external source to balance 
the energy transferred within the mixture to the colder fluid. The 
cold fluid is hotter than its surroundings, to which it loses 
energy . 

The characteristic times for the three different processes; 

energy gain by the electrons, transfer to the neutrals and dis- 

sapation to the surroundings, are different. The latter process 

is very slow. The characteristic time for energy transfer to the 

neutrals is of the order of the electron-neutral collision time, 

- 8 

typically ^ 10 seconds. Thus, if the energy transfer to the 
electrons is caused to vary at acoustic frequencies, by varying 
or E, the time delay can be neglected and the energy transfer 
to the neutrals can be considered in phase with the energy transfer 
to the electrons. Since the time for energy transfer to the sur- 
roundings is much greater, the temperature change, and therefore 



11 



the pressure, of the neutrals will be essentially in phase with 
the energy transferred in. This means that a change in either E, 



the electric field or N^, the electron density, can generate a 
time dependent pressure component in the neutral gas which can act 
as an acoustic source. 

Amplification of a wave occurs when the wave causes a change 

in the medium which leads to a change in the v\»ave in phase with 

itself, i.e. ^positive feedback” occurs. This can occur in a 

discharge plasma through the interaction of the pressure wave in 

the neutral gas with the electron energy source terra H . The 

e 

neutral gas density changes in phase with the acoustic pressure. 
Since both the mobility and electron density depend on the neutral 
density, +he ff^^edback can be introduced through them even with no 
change in the electric field. 

In the simplest analysis, the electron density varies in phase 
with the neutral density. If n and n are the perturbations in 
the electron and neutral gas densities N and N , then 



n 

e 




n 



This indicates that the energy transfer to the neutral gas varies 
in phase with the acoustic pressure variation. Since this would 
produce a pressure variation in phase with the original pressure 
variation, it constitutes an amplification mechanism. For net 
growth to occur, the acoustic losses due to thermal conductivity 
and viscosity must be less than the acoustic source gain. 
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A more sophisticated analysis indicates that phase differences 
occur between the electron, ion and neutral density variations xvith 
resultant generation of a space charge electric field. The phases 
of all of these quantities influence the phase of the feedback and 
determine whether the plasma effects increase or decrease the net 
attenuation of the acoustic wave. 

A. CHRONOLOGY OF THE EXPERIMENT 

Interest in plasma acoustics here at the Naval Postgraduate 
School arose out of the study of moving striations. Earlier work 
performed here ivas that of Partlow Cl], Carr eta and Moore [2], 

Wilson [3]? Crandc^ll [4] and Melville [5]* Melville's work 
represented a shift away from the study of moving striations. This 
was prompted by a theoretical paper by Ingard in I966 [ 6 ], which 
predicted amplification of an acoustic xvave in a discharge plasma. 

A more comprehensive theory by Ingard and Schulz followed in 

1967 [7]. 

The experimental investigation described in this work began 
in 1968 as a natural continuation of the earlier work performed 
here. The earlier work here and elsewhere had demonstrated that 
sound could be generated by a plasma. This coupled with the 
prediction of effects on the sound as it passed through a discharge 
led the author to the consideration of a combined experiment. In 
the proposed experiment, a discharge xvould be used to create sound, 
and it would be observed before and after it propagated through the 
discharge. With the proper techniques, the experiment could provide 
information about both the creation of sound and Its modification 



in the discharge. 
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A numerical analysis of the Ingar d-Schulz dispersion relation 



contained in [7] was performed as part of this investigation in 
which the variables were neutral density, neutral temperature, 
electron density, electron temperature and angular frequency of 
the acoustic signal. The attenuation was minimized with no con- 
straints placed on the values of the variables except for upper 
and lower limits. The results showed that the minimum attenuation 
occurred at high neutral densities, high electron densities and 

temperatures, low neutral temperatures and low frequencies. The 

. . . 2 

maximum relative amplir icat ion was also investigated. The con- 
clusion drawn from a study of the behavior of these two quantties 
was that the experiment should test the theory in the medium 

observed by Yatsui, et.al. [8,9] indicated the need to design the 
experiment so that the polarity of the discharge could be reversed 
Also, it was felt that thermocouples were needed to determine the 
neutral gas temperature in and out of the discharge. 



2 

*^Relative amplif ic:at ion^^( A^) is the change 
due to the plasma amplification referred to the 
uation. It is defined as 



fl^u 

Pu 



A 



R 



in attenuation 
classical atten- 



where p is the attenuation coefficient, I Signifies ionized and u 
signifies unionized. A negative value means that an amplification 
is operative in the plasma case. 
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The experiment was designed to use two mobility-limited 
thermionic diodes as microphones to observe the signal at two 
points with the signal being produced by a remote source discharge. 
Such an arrangement would permit a direct determination of the 
attenuation. The theory of the diode microphones was developed 
beyond the earlier elementary . form so that an extensive cali- 
bration was not needed to use it as an analytical tool. Technical 
problems were solved so that the tube could be constructed with 
a smooth surface, free from side arms, bumps and bubbles which 
affect an acoustic signal, in spite of the numerous tube elements 
contained. A circuit was designed to control the discharge cur- 
rent to produce the sound in the desired form. 

Pulsed currents xvcrc used to produce the sound v/ith the source 
discharge. The current and sound waveforms were Courier trans- 
formed to determine their spectra. By a process of comparison, 
extrapolation and comparison, the frequency dependence of the 
acoustic attenuation and the source transfer function were obtained. 
The measured values of the experimental parameters were used to 
determine the predicted values of attenuation from the Ingard- 
Schulz equations. These equations predicted that the plasma would 
not affect the attenuation under the conditions used. The experi- 
mental value for the attenuation co-efficient agreed with the pre- 
dicted value within experimental accuracy although the results 
indicated the possibility of a plasma effect. 
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The results of Yatsui et.al. [8,9] were then reevaluated. 

It was noted that while the relative change in the attenuation 
coefficient was small, since the attenuation was very large, the 
amplification effect was large. In view of these results, it 
was recognized that the quantity to optimize is the difference 
in attenuation co-efficients between the ionized and un-ionized 
cases . 

Earlier a theoretical study of the frequency characteristics 
of a discharge sound source was undertaken by Schulz [ lo] . He 
concluded that the energy transfer from the electrons was the 
dominant source mechanism. The predicted frequency dependence 
is sin ( V n a/c ) /( V TT a/c ) . When the source transfer function was 
exp^^r im^^'nt A 1 I3/ determined, it was found to be consistent with 
theory at high frequencies but not at low frequencies. The low 
frequency effect was attributed to diffusion effects. 

In the following two sections of this chapter, earlier experi- 
ments relating to the creation and m’odif icat ion of sound are 
discussed. Following, Chapter III discusses the theoretical treat- 
ment of sound creation and modification. This order of presentation 
was chosen because it was felt to offer a better insight into the 
aspects which the theory should treat. 

B. PREVIOUS EXPERIMENTAL WORK ON ACOUSTIC INTERACTIONS 

Before considering specific papers of direct interest to this 
experiment, a quick general review of previous experimental investi- 
gations of the acoustic interaction is in order. Interest in the 
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interaction between charged particles and neutral acoustic waves 
is recent. A relationship between moving striations and acoustic 
waves was suggested by Watanabe and Oleson [ll] in 1955. 

Wojaczeck [ 12] reported observations of a sound interaction in 
i 960 during experiments in Xenon. He noted on effect in his main 
discharge which occurred when an auxiliary discharge was pulsed, 
and he identified the effect as due to sound waves generated by 
the pulsed discharge. 

Since that time, many workers have reported observations of 
sound waves generated by discharges. The source and type of 
disturbance in the discharges, which were observed to produce 
sound, have been sufficiently varied to conclude that it is a 
general phenomenon. Sound has been observed [ 1 - 3 ] whose source 
was interpreted to be the current fluctuations associated with 
moving striations. Sound production by artificially introduced 
current fluctuations has been reported C4,13“l5] for a variety of 
conditions. Modulation of the electron temperature in a discharge 
by an R.F. field was observed by Fitaire and Mantei [16] to produce 
sound. Multiple pulse emission from an HCN laser, reported by 
Turner, et.el., [l7], was interpreted by Born C18] as due to sound 
production by the discharge. 

The variety of techniques used to study the sound produced by 
pulsed discharges is almost as interesting as the observations 
reported [ 4 , 5 , 12 , 19-25] . The pulses ranged in duration from one to 
fifty microseconds [ 23 , l 8 ], in pressures ranging from 5 ^ 10 torr[l9] 
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to 20 torr [l 8 ]. In some experiments, the sound was detected by 
observing its effect on some property of the discharge which 
produced the sound. An example of this is the work reported by 
Berlande, et.al. [24], in which they observed perturbations in 
the light of the afterglow of the plasma which produced the 
sound. Nygaard [21] utilized, the light perturbation effect by 
using one discharge to produce the sound while he observed the 
afterglow of a second discharge. Goldstein et.al. [23] used a 
mobility-limited thermionic diode as a microphone. These and 
other techniques were used in some cases because of the problems 
associated with conventional acoustic detection techniques under 
the conditions involved. In other cases, they were used because 
the sound and ats detection was an une'cpoctcd result in an 
exper iment . 

Only a few papers have been published concerning the effect 
of a plasma on sound propagating through it [5,8,9,26], as 
opposed to the effect on the discharge. Because of their relevance 
and sparsity, these papers will be discussed in detail in the next 
section. 

C. SPECIFIC EXPERIMENTS OF INTEREST 

The two aspects of this experiment are the creation of sound 
and its modification in a gas discharge. Earlier experiments have 
examined both of these and as the experimental results have been 
reported, theories have been developed to explain the experimental 
results. The theoretical aspects of the experiment will be dis- 
cussed in the next chapter . 
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1 . 



Sound Production 



A review of the existing experimental evidence concerning 
sound production can point out features which must be explained 
by any theoretical model of the process. 

Wilson [3] determined from phase measurements of the sound 
produced by moving striations, that the only sound source position 
consistent through the experiment was at or near the cathode. By 
comparing the current and sound waveforms, he concluded that the 
sound was produced by fluctuations in heating processes occurring 
near the cathode due to current fluctuations. He suggested that 
the mechanism must involve an interaction between the fluctuations 
in ion density and temperature and the neutrals. 

Ny gaar u and Meltz l 21] axe mux e specific a^a to the location 
of the source. They state that the acoustic disturbance from their 
H.V. pulse discharge originates in the negative glow region. They 
postulate a high pressure driver region near the cathode maintained 
by electron heating and momentum transfer from the ions. 

Hayess [ 20 ] also interpreted the sound production as arising 
from localized heating of the neutral gas. However his experiment 
brings in a different aspect. He used a high voltage pulse applied 
to an electrode near a discharge already running, so that the other 
electrode of the circuit had a much greater charge carrier con- 
centration than the exciting electrode. He found that the sound 
source was always at the exciting electrode, regardless of the 
polarity of the voltage pulse. By varying the rise time of the 
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voltage pulse, Hayess established that it had particular importance 
in the sound excitation process. Pulses with rise times of 
5 X 10 ^ V/sec produced sound, those with rise times of 1.4 x lO^V/sec 
did not. In addition, sine waves with an average rise time of 

7 

2 X 10 V/sec did not lead to sound production. Hayess also 
established that the presence of charge carriers at the exciting 
electrode reduced the sound produced. From consideration of these 
factors, Hayess concluded that energy gained in the high field by 
the charge carriers is transferred to the neutrals. He further 
concluded that the ions were responsible for the inertial phenomena, 
i.e. dependence on rise time. 

In contrast to the above experiments, the observation of 
soijnr] production by modulation of the electron temperature in a 
discharge, reported by Fitaire and Mantei Cl6], may be deemed as 
experimental confirmation of the roll of the electron temperature 
in the acoustic source term of Ingard and Schulzes equations [?] • 

In interpreting these experimental results, it must be borne in 
mind that there is a distinction between the energy of electrons 
in a low field region, such as the positive column and a high 
field region such as the cathode fall. In the latter, the motion 
is highly directional and can not be considered thermal. 

Upon considering the above experimental results, the 
following conclusions may be drawn. There are two mechanisms 
through which sound is produced in a discharge. One mechanism 
occurs near the cathode and probably involves ions in a high field, 
the cathode fall. Nygard*s contention that the source is the 



20 




t 



negative glow region is in contradiction with the results of 
Hay ess concerning the high field. The other mechanism occurs in 
the positive column and involves the electron temperature. The 
former mechanism appears to dominate the production process. 

Further comment shall be reserved for a discussion of sound pro- 
duction in Chapter II. 

2 . Sound Modification 

The experimental evidence concerning the effect of the 
plasma on sound is not so conclusive. Wojaczek [l2] observed the 
perturbation, produced by the sound, in the light emitted from a 
Xenon gas discharge. On the basis of his calculations, he con- 
cluded that the attenuation was much greater than expected. He 
axtribured this to poijsiblt; interaction Vvith the discharge - 

Later Hayess [ 19 ] pointed out that the attenuation 
Wojaczek predicted was in error because the dominant loss 
mechanism, wall effects, was neglected. Hayess recomputed the 
predicted values. With the wall losses taken into consideration, 
the experimental values agreed with the predicted values within 
the experimental uncertainty. Wojaczek *s experiment was. conducted 
at pressure from .07 torr to 1.5 torr with currents from 1.5 to 
0 amper es . 

Working in Argon, Yatsui et.al. [8,9] studied the effects 
of a discharge on ultrasonic sound. They used barium titanate 
transducers with resonant frequencies of 50 kHz, 100 kHz and 
200 kHz for tr ansmitt€?r s and receivers. The signal received with 
no discharge, as a function of transducer input voltage, was the 
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reference. By dividing this into the results obtained with the 
discharge on, a ^spatial growth rate'* was obtained. A definite 
change in the growth rate vvas observed. It was found to be de- 
pendent on the current and pressure of the discharge and on the 
polarity with respect to the propagation direction. The polarity 
effect was so strong that under conditions such that while propagation 
in the direction of electron drift motion enhanced the signal by 
50 %; no signal v/as received when propagated contra to the electron 
drift. The variation of the growth rate with pressure and current 
was not raonotonic, as the signal exhibited a maximum enhancement 
at intermediate values of the range investigated. The conclusion 
drawn by the authors was that strong damping occurred for contra- 
propagating signals under all conditions. For parallel propagation, 
there is a critical pressure p^ above which attenuation occurs and 
below which amplification occurs. They attributed the directional 
effect to coupling between the sound and the drift motion of the 
electrons, a coupling v;hich is strongest when the sound velocity 
is nearly equal to the electron drift velocity. 

Their experiment provides experimental support to the con- 
cept of energy transfer from the electric field to the neutrals 
via the electrons. It should be pointed out that the acoustic 

frequency is of the order of the collision frequency at the low end 

“2 

of their pressure range, 10 to .5 torr . Indeed it was at this 
end that the maximum effect was observed. 

Macon [26] used a Helium discharge to investigate the 
effect of the discharge on sound. He employed resonant cavity 
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techniques with frequencies below 2000 Hz. The sound source was 
calibrated to provide constant amplitude deflections over a range 
of frequencies. A rigidly mounted microphone received the signal 
at the other end of the tube. By comparing the ratio of pressures 
at adjacent pressure maxima and minima, and assuming the net loss 
was dominated by boundary effects, the attenuation constant was 
determined. The value for the "discharge off" condition was com- 
pared with the "discharge on" condition. Measurements for the 
discharge off condition were taken one-second after turning the 
discharge off in order to have comparable temperature. The tem- 
peratures, as determined from sound velocity measurements, dropped 
from 450 ^K with the discharge on to 360^K ivith the discharge off, 
with one-‘^econd. The results indicated a decrease in alleijuation 
rate of 3^%; a value too large to be accounted for by the pressure 
and temperature differences. Although the experimental results for 
attenuation and amplification rates were both in poor agreement 
with theoretical predictions, Macon concluded that amplification by 
the electrons was definitely occurring. 

These experiments indicate that there are amplification and 
attenuation effects. Wojaczek's experiment does not hint of any 
directional effect. Macon did not look for any. The directional 
effect observed by Yatsui is explainable by models for other fields, 
namely electronics and solid state. The Ingar d-Schulz model to be 
discussed in the next section does not take directional effects 
into consideration . 
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1 1 . THEORY 



Physical theories go hand in hand with experiments. Theories 
are developed to account for experimental results. Experiments 
ate conducted to evaluate theories by testing their predictions. 
In the previous chapter, the experimental evidence regarding pro- 
duction and modification of sound by a discharge plasma was re- 
viewed. In this chapter, the theories will be examined. 



A. SOUND PRODUCTION BY DISCIEXRGE PLASMAS 

After several investigators reported observation of sound 
produced by electrically modulated discharges, Ingard [6] performed 
a theoretical analysis of the problem. He derived the wave equation 
for the pressure in the neutral gas in which the interactions with 
the electrons and ions were expressed in terms of sources in the 
wave equation. Neglecting the fluctuations in neutral particle 
density resulting from the imbalance of ionization and recom- 
bination rates, Ingard obtained the following wave equation: 



1 

2 

c 




V p 



y-1 bH 
2 bt 

C 



V-F . 



( 2 . 1 ) 



In this equation, p is the acoustic pressure, defined as the time- 
dependent pressure variation about the ambient pressure; c, the 
speed of sound; and y is the ratio of specific heats. ”H” is the 
energy source term, representing the rate of energy transfer to the 
neutral gas from the charged particles per unit volume. The rate 
of transfer of momentum to the neutral gas per unit volume is re- 
presented by the momentum source terra F. 
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Ingard considered the effect of momentum transfer represented 
by V*F as negligible for two reasons. First the electron and ion 
contributions tend to cancel each other. Second, only the spatial 
variation of F leads to sound, and for most cases, this is very 
small. For the energy source term, he wrote 

N N (o) (2.2) 

where v^ is the average electron velocity and (or) is the average 
collision cross section for elastic scattering, Ingard considered 
the ion contribution negligible due to its much lower mean energy, 
Ingard pointed out that since the electrons were in an equilibrium 
state, the energy lost to the neutrals, represented by H, ivas the 
same as the energy gained from the field; Therefore 




H = N^eE(fi^E) (2,3) 

where u is the electron mobility and E the electric field, 
e 

Schulz [lO] used equation (2.1) ’as a basis to develop an 
expression for the sound field produced by current-modulating a 
discharge operating between two planar electrodes. He considered 
the effects of the two source terms separately in calculating the 
pressure response. For both, a one dimensional problem was 
assumed and it was assumed that the time dependence was given by 
exp( -iCjOt ) . Under these assumptions, the wave equation inside the 
source region becomes 
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If H is assumed to be spatially constant within the region, 



then the general form of the solution inside of the region is 
_ ^ -ikz ^ ikz . y-1 

P. = Ae + Be + i H. (2.5) 

in CO ^ ' 



With the additional assumption of symmetry about the origin, 
chosen to lie between the two electrodes which are separated by 
a distance a, Schulz’s expression for the pressure field outside 
the source region for the energy source is obtained. 



Pe = C^) 



( 2 . 6 ) 



Schulz assumed that the momentum transfer acoustic source 
arose from an unbalanced fluctuation between the forces exerted 
by the ions and the electrons on the neutral gas. He stated that 
this force per unit volume was given by 



F = N.CO. 



qE ,, CO qE 

- — r— - N en^ 

1 in CO. -ICO e 



in 



CO -ICO 
en 



(2.7) 



with the assumptions CO. >^C0 and CO >^C0. , this simplifies to 

^ in en in 

F iN^qE(CO/OJ^^) . (2.8) 

If the electric field E is assumed to be constant between the 

— > 

two electrodes and zero outside, then V*F can be expressed as 

V-F = F [6(z+|) - 6(z-|)] 

where 6(z) is the Dirac delta function. This leacls to a pressure 
field of: 
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Schulz concluded that the momentum source term would be negligible. 
He compared the two sound fields and found that the ratio was of 

^in A 

the order of — — ) in favor of the energy transfer process. 

" V"/ 

Schultz's results predict that the energy transfer mechanism 
dominates. This result is in agreement with Ingard's earlier [6] 

conclusion. Schulz's expression for the sound field predicts 

. , . . , ^ ^ , . , . , s in X , an 

that the frequency Dependence is like where x = V — . 

X c 

Heretofore this frequency dependence has not been experimentally 
investigated. However there is experimental evidence which refutes 
Schulz's results and contradicts Ingard's assumptions concerning 
the energy source term. 

As mentioned earlier^ the experimental evidence suggests that 
the principal source term lies near the cathode and involves the 
field. Schulz's analysis implied that the positive column was 
the source region. The frequency dependence also shows 

that the sound produced by high frequencies will be small. The 
effective frequencies of the sound pulses mentioned earlier have 
been in the range of 20 kHz to 1 f4Hz , all of which are high 
acoustic frequencies. Perhaps the comparison with pulsed dis- 
charges is not a valid one because of the difference between the 
experiments and the conditions assumed for the derivation. The 
derivation presumes small perturbations yet the current changes 
in a pulsed breakdown discharge can not be small perturbations. 

Since sound production appears to occur in the vicinity of the 
cathode and appears to be associated with high fields, a possible 
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source term compatible with experimental results would be a field- 
dependent energy source term. A field-dependent source is one in 
which the energy for the acoustic wave comes from the field. The 
dependence on this term enters through the energy equation in 
contrast to Schulz’s treatment, in which the field entered through 
the momentum equation. 

When the processes occurring in the region of the cathode are 
considered, a mechanism for the production of sound suggests 
itself to the author. The contribution of ions to the energy 
transfer is usually neglected. Due to their higher mobility the 
electrons transfer energy from the field to the neutral particles 
much more readily. It should be recognized that due to their 
common masses, the ions and neutrals vjill be at es‘=^ent ^ ^ y 
same temperature. The unity mass ratio means that during a 
collision, an ion will transfer a large fraction of its energy. 

In the cathode region, this factor becomes important. There the 
ions are the majority carriers of the discharge current. Their 
density is higher and they accelerate toward the cathode under the 
high field of the cathode fall. The increased energy will be 
transferred to the neutrals through collisions. There will also 
be a net momentum transfer since the ions acquire momentum from 
the field and transfer it also. 

A rough approximation for the strength of such a source term 

can easily be developed. The average energy gained from the field 

between collisions is eE^ . where E is the electric field strength 

in 

and -C is the mean free path of the ions. Since approximately 

half of the ions energ\/ would be transferred in a collision, the 
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energy transfer rate per unit volume is 






H. 

1 




CO. 

in 



N. . 

1 



( 2 . 10 ) 



This approximation is low because it uses the average energy that 
would be transferred if the ions were stopped after a collision. 
It does provide an order of magnitude estimate. If f . is 
approximated by then 

H. <=« eE ”.N. . (2.11) 

1 2 1 1 ' ' 



The 


ratio of this term, 


wher e 


E = B , 
c 


the 


field in 


the 


cathode 


fall to 


the electron source 


term , 


E = E , 
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the 


field in 


the 


pos it ive 


column 


is 















H. 
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%gE v.N. 

N gE (}i E ) 
e pep' 



( 2 . 12 ) 



Now gN M* E is the cuirrent density in the positive column, and 
e e p 

e is approximately the current density in the cathode region. 

Since the two current (densities should be the same, 



H. E 

1 ^ _c 

H 2 E * 
e p 



(2.13) 



The approximations used in arriving at this result are equivalent 
to saying that in the positive column, all of the power delivered 
goes into heating the neutral gas and in the cathode region, part 
goes into heating the neutral gas . The ions accelerate under the 
action of the field in the cathode fall and gain kinetic energy. 
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The difference between the energy taken from the field and the 
kinetic energy gained represents the energy transferred to the 
neutrals through collisions. Since the current density is con- 
stant but the field varies throughout the length of the discharge, 
the energy transfer rate is a function of the field and the 
fraction of energy lost in collisions. 

B. SOUND MODIFICATION BY THE DISCtlARGE 

The theoretical problem of acoustic waves in a weakly ionized 
gas has been studied by several authors other than Ingard [27-333 • 
Without exception, the work of these authors needs modification 
before it can be used to predict the results of an experimental 
study in a laboratory. This is because they neglect boundary loss 
effects; yet at the low pressures of most discharge experiments 
with typical tube diameters of less than, or of the order of, 5 cm, 
boundary losses are dominant. The work of Ingard and Schulz [ 7 ] 
is the most applicable to an experiment. In addition, it is a 
general treatment and is reasonably susceptible to numerical 
evaluation . 

The mechanism of amplification is embodied in the concept of 
energy transfer from the electrons to the neutrals through col- 
lisions. In the absence of an acoustic wave, the energy transferred 
is spatially uniform and merely serves to increase the temperature 
of the neutral gas. In the presence of an acoustic wave, the energy 
transfer becomes time dependent and spatially non-uniform. This 
produces an acoustic v/ave which can be in phase with the initial 
wave, leading to amplification. 



Ingard^s initial analysis of the problem [6] followed the 
experimental observations of sound produced by electrically 
modulated discharges [2,12,13,24]. He treated the problems as 
belonging to two categories, sound produced by current fluctuations, 
or an acoustic wave instability. For the first category, he derived 
equation (2.1). To treat the second category, he placed a term 
representing the electron density perturbation into the energy 
source term of equation (2.1), His results indicated that the 
acoustic signal from a current fluctuation and also the growth 
of an acoustic fluctuation would be small at room temperatures but 
enhanced at cryogenic temperatures. 

Although published in May 1966, the development of reference 6 
was done in early I 965 . In xhe fall of i960, Macon l 2 u] compleLeu 
his investigation. A thesis student of Ingard*s, he attempted to 
conduct his experiment at cryogenic temperatures but had to use 
room temperatures because of noise due to moving striations. He 
concluded that amplification was taking place at temperatures 
greatly above room temperature (450^K). His values of growth and 
attenuation rates were not in agreement with theory but he attri- 
buted this to his experimental arrangement. The more comprehensive 
treatment of the problem by Ingard and Schulz [ 7 ] was published in 

1967. 

The Ingar d- Schulz analysis is based on a three fluid model of 
the ionized gas with both energy and momentum transferred between 
the fluid components. The effects of viscosity and heat conduction 
within each component and v/ith the walls are also taken into 
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B 



consideration. The conditions assumed in making approximations 
may be considered typical for a gloiv discharge: gas pressures of 



The equations forming the basis for the analysis are introduced 
as linearized space-time Fourier transforms of the continuity 
equation (2.14), the momentum transfer equation (2.15), the energy 
balance equation ( 2 . l6 ) and the equation of state (2.17) for each 
fluid component (these equations are presented here for complete- 
ness and all items will not be defined at this time). 



10 12 -3 

1-10 torr, electron densities of 10 - 10 cm , and acoustic 

frequencies 10 ^ CO ^ lO"'^, 



icon . + ikN .V . = 0 
J J D 



(2.U) 




(j,k,f = n,i,e and permutations) 





(2.17) 



- iCO c 0 



ve e 



ik(KT /m )v - / ^ (^?^) ^ 9 

^ e e 2e^ ’ ^ ve e 





(2.18) 



( j=n,i,e) 



The frequency CO is understood to be real and positive with the 
waveform taken as exp ( ikx- icot ) . This means that a sign difference 



between the real and imaginary components of the wave propagation 
constant k corresponds to amplification. Amplification occurs when 
the amplification mechanism dominates and net growth of the wave 



occurs . 

In the momentum transfer equations (2.15), the term /^^. ( k,CO)v^. 
represents the viscous stress, in the j-th fluid component. A 
similar term f^^{kjOO)Q^ accounts for heat conduction in the j-th 
component, where is the temperature perturbation. Poisson ^s 
equation is used to relate the electric field perturbation to the 
density perturbation . 

Ingard explicity assumed that inelastic collisions and ion- 
neutral energy transfer could be ignored. Other assumptions are 
implicit in the analysic because of neglected terms j the d.c. 
electric field is not present in the energy balance equation, and 
ambipolar diffusion to the walls is neglected. 

By algebraic methods, the equations were reduced to a set of 
three homogeneous equations for the component velocity perturbations, 
written here in matrix form. 
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V . 
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L- 



= 0 



(2.19) 



Approximate expressions for these co-efficients were given by 
Ingard and Schulz, in the main body of their paper, as being 
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consistent with the acoustic mode of motion in a gas with a low 
degree of ionization. These equations are: 



A 



D 



G 



1 




iCO. /CO; E ^ 
in 



iCO /CO; H ^ 
en 




Using the co-efficients given in equation (2.20), Ingard and 
Schulz made additional approximations and expressed the dispersion 
relation in the form: 




where 

Cl 

n 



(CO .+ CO ) 6 

' ni ne 





CO 



Cl 

n 






r y ^ 2n 

n ' n 






+ (cor|/2pr^) + (ojic /2/)c r^) 



It was stressed that this expression was valid only if the 
second and third terms in the brackets were small with respect 
to unity . 
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In the above expressions, the adiabatic speed of sound for 



the j-th component, c^. , is defined as: 

% 

c . = (y .KT ./m .) , 

D 3 3 3 

where: V. = c ,/c the ratio of specific heats; K, the Boltzman 
J pj vj ’ 

constant; and m^. are the temperatures and mass of the j-th 

component. The term ^ constant of proportionality 

v/hich is related to the effective collision frequency of a 
particle of the j-th type moving through the fluid component f 
with particles per unit volume. 

The Kirchoff wave attenuation caused by viscosity and heat 
conduction, which was introduced through the terms /_ and /_ 
in a neutral gas is taken into account by use of the character ist ic 



frequency defined by 



CO 

CO 

nn 



1 

cor 




+ 



r 




) . 



( 2 . 22 ) 



The characteristic frequency for ions, is defined in a 

similar manner. The ion and electron plasma frequencies are 

denoted by CO. and CO . The term l/cOT expresses the effect of the 

heating of the neutrals by the electrons. 

In the analysis of their results, Ingard and Schulz point out 

that wave amplification will occur if the terra in the imaginary 
2 

part of k , Eq. (2.21), involving T dominates. They give a simple 
argument and develop a criterion for amplification involving the 
electron-neutral collision frequency 





> 



1 ^ 

r CO 



n en 



(2.23) 
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They also point out that contrary to the result obtained by 
Ingard [6] earlier, this more general treatment indicates that a 
decrease in neutral gas temperature does not always increase the 
amplification. They show that the optimum temperature can be close 
to room temperature, based on the approximate expression Eq. (2.20). 

Solving equation (2.19) for the ion and electron perturbation 
velocities in terms of the neutral velocity and using approximations, 
it is shown that the assumption often used, that v , v. and v are 
all the same in the acoustic mode, is not generally valid. It is 
true only for frequencies below some characteristic frequency. At 
higher frequencies, the ion and electron velocities decrease in 
amplitude and are brought out of phase with the neutrals. 

Maccnc observation of an amp! t cat ion effect at increased 
temperatures does not contradict the Ingar d-Schulz development; 
the observations of Yatsui et.al. do. Since Ingar d-Schulz did 
not consider directional effects, and did not set up the equations 
so as to include any drift motion, their results could not have 
predicted Yatsui *s results. The question arises whether the 
directional effects are present at lower frequencies and higher 
pressures. In the Yatsui experiment, the frequencies were high 
and the pressures low, resulting in conditions v;here the wave- 
lengths were on the order of the mean free paths. 
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Ill . NUMERICAL EVALUATION OF INGARD-SCHULZ EQUATIONS 



The despersion relations developed by Ingard and Schulz predict 
a change in attenuation, possibly even amplification for an acoustic 
wave in a weakly ionized gas. To improve the chances of observing 
what might be a small effect, it was felt that the optimum con- 
ditions must be determined and the experiment conducted under 
similar conditions. '^Condit ions here refers to the values of 
•controlable discharge parameters. 

The purpose of the experiment was to measure a change in the 
attenuation constant. When the terms in Ingard and Schultz's 
expressions were analyzed, the obvious variables were the micro- 
scopic parameters: neutral narticle densitv N . electron density N , 

r ir ^ ' n ' e 

neutral temperature T^ , and electron temperature T^ and the angular 
frequency of the acoustic wave CO. The problem of determining the 
optimum conditions was thus the mathematical pr obLem of minimizing 
a general function of five variable? with the variables constrained 
to a certain region, the region of experimental feasibility. This 
is best done with a computer . 

Since as presented, the Ingar d~Schulz equations were not 
suited to such an analysis, the process was one of analysis and 
synthesis. Analysis was needed to express the terms as relations 
involving only the above variables and physical constants of the 
gas such as cross section or molecular diameter. These terms were 
then evaluated numerically and combined with each other to compute 
the predicted value of the attenuation constant for the input 
values of the variables. 
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Using a computer, the equations were evaluated far a wide 
range of values of the input variables. This provided a more 
thorough evaluation of the Ingar d-Schulz equations than they had 
made. The results showed that their approximate expression, 
equation (2.21) predicted net amplification over a large region of 
variable space although the more complete exj^ressions never did, 
even when Ingard^s criterion for the validity of the approximation 
was satisfied. 

The process of analysis is described in Section IIIA following. 
The procedure used in evaluating the terms and determing the 
optimum values for the variables is described in section IIIB. The 
results of the evaluation are presented and discussed in section 
lilC. - 

It should be pointed out that this section discusses the 
evaluation of the Ingar d-Schulz equations before the experiment 
was conducted. Afterwards, the program was used to compare the 
predicted values of the attenuation constant with the observed 
values by substituting known or estimated values for the variables. 

A. ANALYSIS 

1 • Obtaining the Attenuation Factor 

The system of homogenous equations (2.19) has a non-trivial 
solution only if the co-efficient matrix is singular, i.e. its 
determinant is zero [341. We can therefore write: 

11 co-efficientsll = A(E1-HF )+B(GF-DI)+C{ DH-GE) = 0 (3-1) 
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We define the complex quantity n = ck/co. It can be considered 

2 

an acoustic index of refraction. Since in Eq. (2.20) n appears 
only in the co-efficient A, we may solve for it by defining as: 



A^= A+n 



^ 1 +i 



CO 



CO 



nn 



60 .+60 > 
ni ne 

60 



and substituting this into equation (3*1). After a rearrangement 
2 

of terms, n is found to be: 



2 . CfDH-EG) B(GF-DT) 

n = A + ’ — ^ — — , 

1 EI-HP' EI-HF 



(3.2) 



2 

Now n can be written in the form 
2 

n = \+i)Ji real) 

and the desired result is the wave propagations constant k in the 
fr om: 

k = OJ+ip, real) (3-3) 



so that the problem reduces to solving for the square root of a 

complex number. In particular, since O' is defined to be positive, 

the root of interest is the one for which p has the same sign as . 

2 

This means that if n is written in polar form: 

2 , 2 1 iO 



n = I n I e , 

the solution is the root at Q/2, i.e. the principal root. The 
attenuation factor, p, the ultimately desired result, is given by 



P = f lra(np 

2 

where n^ is the principal root of n . 



(3.4) 
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2 . 



Acoustic Losses 



An acoustic wave is damped through two physically different 



mechanisms, viscosity and thermal conductivity. The damping can 
also be distinguished by the location where the losses occur, i.e. 
at the boundaries or in the bulk of the gas. The latter charac- 
terization has the advantage that it divides the losses into two 
terms of different frequency dependence. When the losses are high, 
the two effects are not simply additive but instead are coupled 
together in the differential equations for an acoustic wave. As 
this would only be the case for viscous fluids or for tubes of 
small diameter^ neither of which are treated here, we shall con- 
sider the effects to be additive. Since the mean free path and 
therefore the boundary layer is inversely proportional to density, 
at sufficiently low pressures the assumption of additive effects 
may have to be re-examined. 

A plane wave traveling in a duct attenuates expont ent ia lly 
with distance with attenuation co-efficient ^^a^\ With small losses, 
a is the sum of the bulk attenuation co-efficient a^ and the wall 
attenuation co-efficient a . For a circular pipe, these co- 

CO xr I- 7 

efficients are given by Morse and Ingard [35jp989] as: 



^CO 2r 



CO/c 



Nj 




(3.5) 
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P 



where K, the thermal conductivity, is given by: 



3 



and 



= LBs. i 

\ 

The mean free path i is 



^ z: — t — 3 L 3 -the co-efficient of viscosity, 



I = 



f2 ttd'^N 



(3.6) 



In these expressions, p is the mass density, c the adiabatic speed 
of sound and D the atonic diameter. When these expressions are 
used for a three component fluid, c will be written c^; however 
there is no difference between the two, c ^ c. 

Sincp tho energy of a wave is given by the square of its 
amplitude, the energy of a wave attenuates with co-efficient 2a. 
V^ith a velocity of c^, the characteristic time for energy dis- 
sipation by the neutral component, T^, is related to the atten- 
uation co-efficient a by: 



y 2a c • 
n 

Therefore, for the small damping approximation 
in a circular pipe, the characteristic time for 
is given by^: 



(3.7) 

of a plane wave 
ener gy dissipation 
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pc CO 
^ p 



(3.8) 



3 . . . . 

The expression does not exactly agree with the expression in 

Ref. 7 but it is consistent with the derivation there and agrees 
with the same expression in Morse and Ingard [351. The difference 
may be due to an error in typesetting. 
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Substituting for K and 77 , using y = 5/3 (a good approximation for 
a noble gas), the above expression simplifies to: 



^ = \F^ 



2 i 5 

CO — + 

C 3r 

n 



J2.4 



cof c 



(3.9) 



This term enters the dispersion relation as CO/cO defined by 

^ nn 

equation (2.22), 

CO ^ 1 

CO cor ' 

nn n 

3 . Zero Ionization Limit 

In order to be consistent, the attenuation predicted by 
Ingard's equations, for the limiting case of zero ionization, 
should be comparable with previous results for acoustic atten- 
uation. For the standard ot comparison we snail use tne 
attenuation co-efficient ^^a'^ which was used as an input in these 
equations . 

a = a + a 

B 






2 

2 ^ 



2.4 CO 



(3.10) 



n 



For zero ionization, equations (2,21) and (3.2) reduce 



to the same expression for n , 



2 , .1 
" ‘ ^ * ^—T ■ 



(3.11) 



Now, both equations (2.21) and (3*2) were based on 
approximations. When the expressions for the co-efficients, 

Eq. (2.20), were presented, Ingard and Schulz had already made 
several approximations. Going back to the full expressions for 
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the co-efficients, which are contained in the appendix of Ingard 
and Schulzes paper, an examination shows that when the ionization 
goes to zero, only one independent co-efficient remains, given 
by : 



A = 1 
o 



c k\2 



1- 



y-1 ^ 



y l + i/' /CO 



f 



+ 1 



In 

60 



Even though there is only one co-efficient, it must equal 

zero for the equations to have a non-trivial solution for V . 

2 

.Setting A^ equal to zero, and solving for n = — - — , we obtain 



CO 
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n = 



1 + i f, /60 

In 



1 - 



y-1 ^ 2n^'^ 



( 3 . 12 ) 



where /. and /' are the viscous and thermal conduction com- 

ponents of the acoustic losses given earlier in section IIA. 

They are hidden in the expression for the attenuation co-efficient, 

The attenuation time T can be expressed in terms of these com- 

n ^ 



ponents as 



-L. = ;• . Zzi / 

r 'in y ' 2n 
n 



Using the same expressions and approximations as before 
they can be written as; 

^ In 



^ I + 7\l\l2,4 co.«c 



2n 
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3 C 3 
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since there were fewer approximations involved in 
Eq. (3*12), the value of p obtained is a stronge): comparison test 
then the use of equation (3.11)* Using both equations (3*11) ^md 
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(3.12), B was calculated for several different values of N , T 
and CO. The results were compared with the value of ”a^^ calculated 
for the same input values using equation (3*10). In all cases, the 
results agreed within less than 1%. It is therefore concluded that 
in the limit of zero ionization, the Ingar d-Schulz dispersion 
relations predict the correct attenuation. 

4 . Neutral Heating by Electrons 

The effect of energy transfer to the neutrals by the 
electrons is expressed through the term 1/oOT. Ingard and Schulz 
define T by the relation: 

i = (y-l) . (3.13) 

C 



The rate ^f energy transfer per unit mass, H, is defined by: 



H = 
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(3.14) 



where O’ , the average cross section for energy transfer, is 
en 



def ined as : 



,3 p/kt 
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(3.15) 



The integral over cos 9 is the total cross section 
(velocity dependent) for momentum transfer, as defined by 

present [ 36 ]. It is given by: 



This is the correct expression for O , vice the expression 
given in [7] 21 s confirmed by private communication with M. Schulz 
(due to a typographical error). 
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(3.16) 



Therefore: 
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(3.17) 



A Maxwellian distribution of electron velocities is a 
good approximation to the actual distribution in a gas discharge 
only at extremely low values of E/p C39J. It is recognized that 
E/p is not low in the usual case as evidenced by figures 125 and 
126 of von Engle [4o]. However the Maxwellian distribution is 
computationally convenient to use and is valid as a first approxi- 
mation to the real distribution function. It is felt that the 
improvemerit in a.c curacy resulting from the use ot a better 
approximation such as the Druyvesteyn distributicn does not justify 
the increased computational difficulties. The value obtained for 

a under this approximation is smaller than the actual value 
en 

because the approximation indicates a greater number of high 
energy electrons than actually exist and the cross section de- 
creases at higher electron energies. Substituting the expression 
for the Maxwellian distribution into Eq. (3.17) gives: 
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The cross section for momentum transfer is available in 
terms of electron energy expressed in electron volts. By defining 
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Et = KT^/e, 
effecting a 
simple form 



the energy equivalent election temperature, and 
change of variable, expressed in the 
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where R = -rr;::;: — is the variable of integration by which the 
2K.J 

e 



average cross section for energy transfer 
the momentum transfer cross section data, 
5* Collisional Losses 

If we refer back to equation (2.2) 
involving the constants of proportionality 
collision integral for momentum transfer, 
can be written: 



can be evaluated from 



, we see that the terms 

CO., arise from the 
3k 

The collision integral 
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bf 



m . ( V .-u ) 

3 3 3k' 



/ dv .= p .to ( V ,-V . )-l-p .CO ., ( 
/coll ^ ^ ^ 



V - V . ) 



In order to give a numerical value to i't must be expressible 

in terms of the variables and gas properties. To do that, its 
definition must be understood. Since tlie momentum transferred 
from species 1 to 2 through collisions is equal and opposite to 
that transferred from 2 to 1 through collisions, then we can 
equate: 



This shows that the constants CO., are related to CO, . as : 

3k kj 

“k 

CO. = — — CO, . . 

ik N . m . kn 
3 3 



(3.19) 
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Notice that due to the mass ratio factor, CO. 



jk 



and CO, 






cannot both 



be collision frequencies. If in Eq. (3.19) the j-th component is 
taken to be the heavier component and the k-th component the ion 
or electron component, the resulting expressions will be con- 
sistent with the relations given by Ingard and Schulz. In 
addition, when the heavier component is the neutral gas, this is 
consistent with the picture of CO^.^ being the collision frequency 
of a particle of the j-th type moving through the k-th fluid 
component since the neutrals are much more numerous. We can there-’ 

fore develop analytical expressions for CO , CO. and CO .. The 
^ ^ ^ en in ei 

other terms will be related by Eq. (3*19). 



section is essentially the same as for a neutral-neutral collision, 
then it follows that the ion-neutral collision frequency is the 
same as the neutral-neutral collision frequency. The neutral- 
neutral collision frequency is very closely approximated by the 
elementary kinetic theory result for a hard sphere ideal gas. The. 
ion-neutral collision frequency can therefore be expressed by: 



b. Electron-Neutral Collision Frequency 

In order to obtain an expression for the electron- 
neutral collision frequency, we must start with the integral form of 



a . 



Ion-Neutral Collision Frequency 



If it is assumed that the ion-neutral collision cross 




(3.20) 



the expression given by Present [ 36 ] for the collision rate for two 



intermixed species. 
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For electron-neutral collisions, v -v ^ v and the integrals 

' e n' e ^ 

can be separated. The frequency of collision of an electron with a 



neutral particle, CO . is 9^ /N . Substituting also the expression 
for the Maxwellian distribution for f (v) , the collision frequency 
can be written: 

mv 

fm \3/2 f CP ■ 2KT 
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dv 



Finally this can be expressed in a form similar to the expression 
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c. The Electron-Ion Collision Frequency 

The electron-ion collision frequency for momentum 
transfer is given by Tanenbaum [4l] as: 



-12 



CO . = 3.62 X 10 N. T 
er 1 e 



- 3/2 



2n -A. (N. in cm"^) (3.22) 



where -A. = 12 n r'l N is the Coulomb cutoff parameter, 

D e ^ 



and r^^, the Debye length is given by: 

. - f j- 
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Expressions for the plasma frequencies for electrons 



CO 



and for 



ions CO. 5 
1 



CO 
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CO 
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are well known and present no computational difficulty. 

6 . lon-Acoust ic Losses 

Losses caused by viscosity and thermal conductivity in 
the ion fluid component are represented by the term CO/CO^^ which 
• is defined by 

/ y . - 1 

if: = 5 -f- • 

The quantity represents viscous stress in the ion fluid. The 

other term f 2 i’ multiplied by 9^, the ion temperature per- 

turbation, accounts for heat conduction in the ion fluid. Adapiting 

a heuristic approach, assume that CO/co. . has the same form as 

11 

CO/CO^^ given by equation (3*9) • Under such an assumption, it can 
be expressed as: 



CO 

CO. . 
11 



f = 



i 






coi . c . 
1 1 



(3.23) 



where the subscripts indicate that these are ion terms. 

It is not unreasonable to assume that, since their masses 
are essentially the same and their temperatures nearly the same, the 
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ion sound speed is the same as the neutral sound speed c^. 
However the mean free path for ions f must be unlike that for 
neutrals, since one must consider the effect of Coulomb inter- 
actions. Wu [42] expresses the mean free path for Coulomb 
interactions, which have been corrected for Debye screening, as: 



L = \ , (3.24) 

N r /x ) 

e m ^ P m ' 

2 

where r = — is the mean distance of closest approach for 

m 4 ttg KT 
o 

two like charges q. Wu refers to T as the mean plasma temperature. 
Since we are considering ion losses, the temperature to be used in 
the expression for r^ is the ion temperature. However the Debye 
length is determined by the electrons due to their higher tem- 
perature and greater mobility. 

Because the mean free path of ions in a weakly ionized 
gas is limited not by Coulomb interations but through collisions 
with neutrals, this expression alone would give erroneous results. 
The mean free path, X, can be expressed in terms of the mean 
velocity v and collision frequency V as 




In a mixture of gaseous components, the total collision 
frequency of one component is the sum of the partial collision 
frequencies. This leads to the result: 




1 
X . 

le 
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Ignoring l/X . ^ and making the approximation that 
mean free path for ions can be written as: 

i = 1 

i l/l^.+l/i ' 

1 

B, COMPUTATION AND EVALUATION 

1. Defining the Limits of the Variables 

Since the goal of the analysis is a prediction of results 
for a given experimental arrangement, it would have been desirable 
to formulate the problem in terms of filling pressure p^, filling 
temperature tube radius r, frequency v and current- I^. Unfor- 

tunately the conversion of p , T , and I to the variables T , T 

^o’ o o e n 

and N^ is not a trivial problem. For that reason, the independent 
variables chosen were: neutral density , neutral temperature , 
electron density N^, electron temperature T^ and angular frequency 
CO. The limits of these variables were then chosen to correspond 
to the values expected in a feasible experiment. 

A practical range of values for the tube radius is from 
.5 cm to 10 cm. This value -was fixed at 2 cm except when investi- 
gating its effects. 

The lower limit of the angular frequency was fixed at 
100 radians/second because it was experimentally possible and 
numerically convenient. The assumption of plane waves was implicit 
in the formulation of the problem by Ingard and Schulz. From an 
experimental point of view, in order to avoid the effect of higher 
modes propagating in a tube, this sets an upper limit for CO of the 



= X , 

nn 



( 3 . 25 ) 
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order of 2nc/r[35]. However most pulse experiments involve 
higher frequencies and Yatsui, et.al, [ 8 , 9 ] worked with 50 > 100 
and 200 kilocycle continuous and burst sineisoidal signals. The 
limit of validity of the continuum description is when the wave- 
length is of the order of the mean free path. This corresponds to 
5 

about 10 radians/sec. at room temperature and pressures around 
1 torr and was selected as the upper limit for CO. 

l4 -3 

The neutral density limits were set at 3*5 ^ 10 cm 

38 “"3 

and 1 . 5 X 10 cm , corresponding approximately to filling pres- 
sures of .01 torr and 50 torr at room temperature respectively. 

The electron-density limits relative to the neutral 
density were selected arbitrarily, with a logical basis, as 

10 to 10 with absol’^te iTTr^-i-fc of lO^nm ^ to 1 . y 10 ^^ cm 

n n ' 

For currents of milliamperes at low pressures around 1 torr, the 
degree of ionization is about 10 ^ [40,43]* A lower limit two 
orders of magnitude below this is reasonable and a higher limit of 
three orders of magnitude greater would correspond to ampere cur- 
rents. These relative limits were tempered by the current densities 
which they implied for large pressure changes. 

The neutral gas temperature is not truly independent 
but actually depends on the discharge current, pressure, tube 
radius, and external cooling effects. For the purpose of this 
analysis, it is sufficient to assume the temperature can range from 
the boiling point of the gas to the upper limit of temperature 
which can reasonably be obtained in an oven, approximately 700 Kelvin. 
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The electron temperature is definitely not an in- 
dependent variable, being determined by the discharge conditions 
specified through the other variables [4o]. The only analytical 
expression available implies that the electron temperature is a 
function only of the gas and the similaiity parameter pR. The 
expression is of limited validity, falling at very low values of 
pR and high p. Rather than restrict Tg with an artificial con- 
straint, it was felt that treating it as an independent variable 
would yield more valid results. The limits were chosen to cor- 
respond to those normally encountered e:q:)er imentally [40,44], 

4 5 

10^ to 10-^ Kelvin. 

It should be pointed out that permitting these 
variables to move freely within their limits places less con- 
straints on their possible values than actually exist. The re- 
sult may predict physically unrealizable results. 

2 . Computation of Attenuation and Determination of Minimum 

For the computation, the formulas w^er e first placed in the 
order in which their values were needed. They were- then translated 
for computer processing and a program written using the FORTRAN IV 
programming language- The values for the fundamental constants 
ate from [45], data for Argon are from Cook [46] except the atomic 
diameter and scattering cross section for momentum transfer which 
are from McDaniel [39]. 

The scattering cross section data were taken from figure 
4-3’*l(c).of McDaniel, using a transparent overlay grid and dividers 
to establish a table of values of versus electron energy’ in 
electron volts. A subroutine was written to use this table to 
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perform the numerical integration to determine using Bq. ( 3 . 18 ) 

and to determine using Eq. (3*21). In computing ^^^9 

momentum transfer cross section data were used to approximate the 
collision cross section. The momentum transfer cross section is 
from 20% to 50 % smaller than the collision cross section for all 
values of electron energy. This substitution causes the electron 
collision frequency to be low by ^ 25 %. If this was found to be 
critical, a correction could be applied. This procedure was con- 
sidered to be justified by the increased convenience in setting up 
the problem and in reducing the number of complete integrations 
required by a factor of two. 

A set of fixed values: CO = 10^, N = 10^^, N = 10^^, 

n e 

o 4 . • 

T = 3^0 K and T “ 10 K. were used to evaluate all t with 

n e 

an electronic calculator. This served as a cross check on the 
computer results. To check the validity of the expressions used to 
calculate attenua^t ion , a comparison was made of the attenuation 
factor for bulk losses with a value from another source. Using 
Eqn . (3.9) ^nd (3.7)? we can write 

^ ^ 

CO 2cC0 I n /BULK c 



ther ef or e 




Evaluated at standard temperature and pressure, this yields 

- ] .09 X 10 m /sec 

/ 
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; for Argon 



2 

which agrees with the value of p V/ j 1.87 x 10 
in Kinsler and Frey [47l . 

3 • Discussion of Results of Calculations 

At the point used for checkout of the computer program, 
the agreement between the attenuation constant calculated using 
(3*1) Ingar d- Schulz approximation, Eq. (2,21), was 

very good with differences on the order of a few parts per million. 
However the agreement was not so good in other areas. The use of 
Eqn . (2.21) yielded predictions of amplification over a wide range 
of conditions. In contrast, results using equation (3*1) show no 
amplification within the entire range of the variables, but merely 
an increase or decrease in attenuation, with respect to that with- 
out ionization. The amplification predicted by Eq . (2.21) is 

therefore considered erroneous, arising from factors being neglected 
in the approximation. 

The minimum attenuation was found to occur at the limits of 
the allowed ranges of the variables. Tliose variables at their high 
limits were neutral density, electron density and electron tempera- 
ture. The frequency and neutral temperature were at their low 
limits . 

Figure 1 shows the attenuation as a function of angular 
frequency, displayed on a semi-log plot. The curve parameter is 
neutral density. Neutral temperature, electron temperature and 
degree of ionization have been held constant at their optimum values. 
A more revealing display is shown in Figure 2 which is a plot of the 
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Figure 1. Predicted Attenuation Factor in Discharge 

Pl^isma in Dc^pendence on CO for Various Values 

of Neutral Density. Parameters are: X = 8 Sk, 

T = 105k and N /N = 10 ~ 5 . " 
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Figure 2, Relative Change in Attenuation due to Ionization 
in Dependence on CO for Various Values of Neutral 
Density . 

Values calculated using Ingar d-Schulz equations 
with parameters T = 88 k, T = 105k, and 
-5 ’ e 

N /N = 10 . The values of neutral density 

e n 

for the curves are: 

1 8 

1. N = 1.5 X 10 

n 

2. N = 2.4 X 10^^ 

n 

l6 

3. N = 3.9 X 10^ 

n 
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fractional change in attenuation due to ionization. The curve 
parameters are the same as for Figure (1). Positive values cor- 
respond to an increase in attenuation; negative values, a decrease 

The ion perturbation velocity was found to remain the same 
magnitude as the neutrals and in phase. The electron perturbation 
velocity was the same magnitude as the neutrals, and in phase, at 
low frequencies, and decreased to zero and went out of phase at 
high frequencies. The frequency at which the transition began 
shifted to higher frequencies as the density increased. This is 
displayed in Figure 3* 

The results of the computation indicate that net amplifi- 
cation is a false prediction, arising from approximations; however 
attenuation is aflected to a measurable degree. The maximum de- 
crease occurs with high pressure, low temperature, low frequency 
and a high electron temperature, highly ionized discharge. 
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Figure 3* Magnitude of Ratio of Electron to Neutral 
Perturbation Velocity in Dependence on 00 
for Various Valvies of Neutral Density. 



Values computed using Inga r d-Schulz equations 

with parameters: T = 88 k , T = 105k, and 
c- n ' e 

N /N = 10 . The values of neutral density 
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IV. EXPERIMENTAL ARRANGEMENT AND PROCEDURE 



The original concept was to use an auxiliary discharge, with 
a gated sine wave of current to produce gated sine waves of sound, 
the sound amplitude was to be measured before and after passing 
through the main discharge area. From a comparison of the signal 
amplitudes at the two points for various frequencies, the attenu- 
ation factor and its variation with current magnitude and polarity 
would be determined. The analysis of a preliminary experiment 
showed that other factors had to be considered and that a Fourier 
analysis of the signals was needed. Compromising betiveen the need 
for a signal with low frequency components and one with all signi- 
ficant features appearing before the reflection of the signal from 
the tube ends, a one millisecond current pulse was selected xo 
provide the sound pulse. Although a two ampere power supply was 
available for the main discharge, the high temperature of the main 
discharge column limited the experiment to half an ampere. 

A. SOUND SOURCE 

The sound source was a small auxiliary discharge far removed 
from the main discharge and the two diodes bounding the main dis- 
charge. The auxiliary discharge was always on when in operation. 

The quiescent state was a normal glow with a current of less than 
10 ma . The peak current during a pulse was around 150 to 200 ma . 
which corresponded to an abnormal glow. The current in the auxiliary 
discharge was controlled through a transistor load. During pre- 
liminary work, it was ascertained that the sound produced 
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closely related to the current in the auxiliary discharge but not 
to the voltage or to the power dissipated. Therefore only the 
current was controlled and monitored. 

Although a pulse was used for the final data, the current con- 
trol circuitry was designed for sinusoidal currents, either con- 
tinuous or gated. To avoid extinguishing the discharge, a minimum 
current flow was required at all times. This required a bias signal 
equal to the sinusoidal signal. The continuous mode was not opera- 
tional due to rapid deterioration (in times less than or of the 
order of one-second) of the discharge current when run at the high 
currents required for a good acoustic signal. The cause was not 
thoroughly investigated, but was believed due to excessive power 
diss ipa t ion . 

The current control circuitry is shewn (enclosed in dashed 
lines) in Figure 4. It consists essentially of two circuits feeding 
the base side of two transistors operating in parallel. One cir- 
cuit consists of a small d.c. battery with a variable resistor for 
controlling the quiescent state of the discharge. The other cir- 
cuit sums two input signals to provide the a.c. signal. 

Because of their high voltage capability (300 V), and their 
high dynamic resistance. Motorola MJ3202 transistors were used. 
Parallel operation was chosen to give a greater current range to 
the discharge while remaining within the power rating and linear 
range of the transistors. Emitter biasing resistors were used to 
balance the load between the two transistors. Although some 
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negative feedback occurred due to the absence of bypass capacitors, 
it did not have a noticeable effect, 

A General Radio GR1217-C Unit Pulse Generator provided the d.c. 
pedestal for the gated sine wave. The degree of modulation by the 
gated sine wave and the amplitude of the a.c, discharge currents 
were adjustable through variable resistors. For ease of control, 
multiturn potentiometers were used for control of bias, quiescent 
current and a,c. current amplitude. 

The General Radio equipment used required the low side of the 
transistor circuit to operate at ground potential. To minimize 
electrostatic interaction with the diod€‘ microphones, the anode of 
the auxiliary discharge was grounded. This meant that both ter- 
minals of the d.c. power s\ipply were floating above ground. Since 
the floating potential varies greatly during the current modulation 
capacitive coupling of the terminals to ground was found to distort 
the discharge current waveform. The po\^er supply used (Universal 
Electronics Model G3050M) was chosen because of its low coupling. 
However its maximum voltage was 3 OO volts. This did not permit 
using the full current capability of the circuit at low pressures. 
Consequently operation was limited to pressures above 5 torr for 
good acoustic signals, although use of another power supply would 
have extended the operating range. Figure 5(a) shows a typical 
discharge current pulse waveform with a gated sine wave. 

The electrodes for the auxiliary discharge were two tantalum 
cylinders’ ( 1^4 cm long by 2 cm dia . ) v;ith faces separated by 3 % cm. 
They were mounted axially in the tube, supported by .070 inch 

diameter KOVAR leads. The leads were beaded with glass for electr • 
cal insulation. 63 
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Figure 5. Sound source current waveform and resultant 
sound. Ambient pressure is 8.7 torr . Gas 
is Argon. Frequency of sine wave being gated 
is 3000 Hz. Time scale is .5 ms/cm. 

(a) Sound source current. Vertical scale is 
100 ma/cm. 

(b) Sound produced by discharge current. 
Vertical scale is uncalibrated. The second 
pulse of sound is the reflection from the 
near end of the tube. 
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B. SOUND DETECTION 



Mobility-limited thermionic diodes were used to detect the 
sound waves. They were mounted axially, one on each side of the 
discharge region. Electrostatic shielding was provided by a 
cylinder of heavy, 2.3 mil., copper foil which was wrapped around 
the tube at the location of the diode. The shield extended several 
inches on both sides of the diode and was connected to ground. 

Since the cables and connectors were shielded, the only part of the 
circuit not shielded was a small portion of feedthroughs where 
they penetrated the glass. Also the diode was not shielded on the 
ends. For a more complete discussion of the diode and its sus- 
ceptibility to noise, see Appendix A. Figure 5(b) shows a typical 
diode response signal. 

For the experiment, filament power for both diodes was provided 
by the same power supply, a Power Designs 36 IOOR regulated power 
supply operating in the constant voltage inode. Figure 6 shows the 
circuit employed for power and signal detection. 

C. TUBE DESIGN AND CONSTRUCTION 

There are two features of the tube used in the experiment which 
are worthy of note. One is the method of mounting the elements. 

The other is the length of the tube and the placement of the elements 
therein. 

The tube elements were mounted on feedthroughs which did not 
distort the tube walls significantly. The feedthroughs were 
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Figure 6. Sound Detection Circuit. 
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sealed with a special glassblowing technique which is discussed 
in the appendix. The result was a tube which could be treated as 
an acoustic waveguide for frequencies much higher than those used. 

To reduce the effects of reflections and electrostatic inter- 
actions, the elements were physically separated so that the acousti 
time of flight separated the unwanted, from the wanted, signals. 
This required a half meter buffer region at the tube ends and a one 
meter region between the sound source and the nearest diode. 

In addition, it was discovered during preliminary experiments, 
that there were positively and negatively charged particles ex- 
tending into the supposedly neutral region far beyond where they 
were expected. A series of experiments established that they were 
real and that both typos were present. The buffer region reduced 
the effect but screening electrodes were used for further reductiorj 
Since with cold filaments, any diode signal was attributable to 
charged particles rather than sound, the screnning effectiveness 
was studied with cold filaments, in which case, the diodes acted 
like double Langmuir probes. The cold filament effect v^’as strong 
in Neon but barely noticeable in Argon. Although interesting in 
itself because these signals were detected at distances from the 
discharge of 30 , 40 and 50 cm in Neon at pressures ^ 10 torr , the 
phenomenon was studied only enough to eliminate its interference. 

A schematic drawing of the tube appears in Figure 7 £^nd the lo- 
cations of the tube elements appear in Table 1. 

Copper -ccnstantan thermocouples, of 5 rail wire, were used to 
measure the temperature on the tube axis at three points. For 
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Figure 7. vSchematic Drawing of Discharge Tube. 

Identification of Symbols and Locations 



are Shown in Table I. 
TABLE I 

LOCATION OF TUBE ELEMENTS 




ITEM 
End ] t 


SYMBOL 


LOCATION 

-149.7 


Anode 1 


Aj 


. 1 r\o 1, 
^ • -1 


Anode 2 


A2 


-98.7 


Thermocouple 1 


TCI 


-50.3 


Diode 1 


D 1 


0 


Screen 1 


SI 


26.3 


Cathode/Anode 1 


Cl 


37.3 


Thermocouple 2 


TC2 


102 


Cathode/Anode 2 


C2 


122.2 


Screen 2 


S2 


133.4 


Thermocouple 3 


TC3 


139.3 


Diode 2 
End 2 


D2 


152.5 
203 + 



* Location given is location of center of element 

T Tube ends taper from 4 cm to 1 cm diameter over 5 cm. 
Distance refers to beginning of taper. 

t For calculating tube volume, the effective length is 3^3 cm. 
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support and insulation, the leads were first sealed in fine twin 
capillary glass tubing with ULTEK VACSEAL epoxy. The wire to 
extend into the tube was left exposed. During construction of 
the main tube, a small piece of tubing, large enough to hold the 
capillary tubing, was attached. The capillary tube was sealed 
inside of this tube with the epoxy. The seal was good and the 
thermocouples worked well. Such a method is not recommended for 
the main discharge area. This is because the temperatures en- 
countered in the vicinity of the discharge may exceed the maximum 
recommended for the epoxy, 150*^C. In this experiment, such heating 
caused contamination of the discharge and limited operation of the 
main discharge to currents below .5 ampere. 

The pumping train used for evacuating the tube consisted of a 
VEECO Model EP2A1 diffusion pump, rated at 85 f/sec, pumping to a 
DUOSEAL forepump. A NRC model 720 ionization gauge control with a 
type 53^“SK ionization gauge was used to monitor the pressure above 
the trap, which was kept filled with liquid nitrogen. A one inch 
line was used between the trap and the high vacuum isolation valve 
(a GRANVILLE PHILLIPS type C value) used to connect the discharge 

tube to the pumping system. The ultimate vacuum obtained was 

- 8 

around 5 x 10 torr as measured by the gauge. 

The gas filling system was connected to the tube via a (CEC 
type 25 ) high vacuum valve. Pressure was measured with a manometer 
filled with Octoil S. A trap located between the filling system 
and the tube was partially filled with liquid nitrogen whenever the 
connection to the tube was open. 
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D. MAIN DISCHARGE 



In order to determine the effect of polarity, the electrodes 

at both ends of tie main discharge needed to serve as cathode or 

anode. Therefore they were of identical design, consisting of a 

filament and a cylindrical anode. The filament was used as the 

hot cathode for the main discharge or as the hot cathode for a 

"Pupps discharge^’ at the anode of the main discharge. It was a 

% inch diameter, 7 turn helix, directly heated tungsten dispenser 

cathode designed for laser work and produced by SPECTROMAT. It is 

rated at an emission current of 5A/cm^ at 1150^C . The anode 

B 

cylindrical electrode which was mounted concentric with the fila- 
ment, was located so that the first two turns of the filament, 
on the discharn<= side, were exposed. This was done so that the 
Pupps discharge would tend to go toward the leading edge of the 
anode, and thus eliminate anode oscillations. A small low 
ripple (20V, 1.5A) D.C. power supply was used for the Pupp^s dis- 
charge power source. 

Filament power for the main discharge cathode w^as supplied by 
a KEPCO KS-8 100m power supply. Filament power for the Pupps dis- 
charge cathode was supplied by a POWER DESIGNS PD 36 IOOR unit. The 
load for the main discharge was two high power variable resistors 
connected in parallel. The main discharge power supply was a 
KEPCO 77 OB, 600 V d.c. supply. 

A reasonably simple, yet effective system was used to elec- 
trically isolate the main discharge system from the remainder of 
the circuitry. The 60 cycle input pov/er for all units was "solated 
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using transformers. All equipment for the main discharge u^s 
mounted on the same rack, which contained no other equipment. The 
chassis ground for all of this equipment was connected to a common 
point which was then connected to ground via a one megohm resistor 
in parallel with a 10 JIf capacitor. The a.c. and d.c. potentials 
across these points were monitored with high impedance devices. 

This precaution was taken to prevent the buildup of a high elec- 
trostatic potential between the two ground points. The a.c. current 
in the discharge was monitored by using a current probe with a 
sensitivity of 1 mv/ma. To prevent a large potential difference 
between the two ground systems, arising from interaction of the 
fields, a potentiometer was used to vary the potential of the main 
electrodes with respect to the secondary ground. Tiiit^ shown 
in Figure 8 as part of the main discharge circuit. 

The discharge was ignited with a Tesla coil in conjunction 
with a copper wire wrapped around the discharge tube between the 
anode and the cathode. This was necessary because of the high 
breakdown voltages required at the pressures used. 

E. ELECTRONIC MEASUREMENTS 

The three waveforms of interest - the sound source current, 
and the two diode signals - together with a trigger signal were 
recorded on a four channel FM tape recorder (Precision Instrument 
Co., pi- 6200) for later analysis. Recording was done at 37-5 ips . 
Record and playback specifications are contained in Appendix C. 

The frequency response did not appear to distort the two diode 
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~ MAIN GROUND 



Figure 8, Main Discharge Circuit 
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signals but it did affect the trigger and discharge current signals 
since these had significant high frequency components. Figure 9(si) 
shows a typical diode signal. Figure 9[b) shows a typical playbacl-: 
of a diode signal. It should be noted that these figures are not 
derived from the same sample. 

The trigger signal was provided by a General Radio GR 1217-C 
Unit Pulse Generator operating at a pulse repetition rate of 
10 Hz. The output signal vjas a one millisecond rectangular pulse 
which was used to trigger the scope and was also recorded. The 
"delayed sync" output was used to trigger another unit pulse 
generator from which a one millisecond pulse was fed to the cur- 
rent control circuit. The discharge current was observed by 
measuring tlie voltage across a 5 ohm precision res is lor. 

The output signals from the diodes were amplified with two 
differential amplifiers, TEKTRONIX RM122, with an input impedance 
of 10 megohms. The amplifiers could provide amplification of 
either lOOX or lOOOx. Roth output signals were recorded on tape. 

The diode output contained very sharp noise spikes which, after 
amplification, adversely affected the recorder operation. To 
eliminate this, a slicing circuit, which clipped all signals whose 
value exceeded I.5 volts, was employed between the amplifier 
output and the tape recorder. This circuit permitted an improvement 
in the signal to noise ratio since the noise could not exceed 
1.5 volts and the signal could be amplified to a level just below 
that . 
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Figure 9. Comparison of Typical Diode 
Signal with Playback of 
Diode Signal . 

(a) Direct display of diode signal. 

(b) Playback of diode signal. Due 
to triggering problems, the 
playback is not a playback of 
(a) . The playback does show 
how little the recording process 
distorted the signal. 
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F. EXPERIMENTAL PROCEDURE 

After preliminary experimental work had established the back- 
ground, the progress of the experiment can be divided into five 
stages: (1) construction of the tube, (2.) preparation of the tube, 
(3) calibration of the diodes (4) data accumulation, and the last 
stage, ( 5 ) analysis of results. The ana.lysis of results will be 
discussed in the next chapter . 

1 . Construction 

Due to its great length, the tube ivas constructed in three 
sections. After construction and before annealing, each section 
was leak tested with a helium mass-spectrometer. The sections 
were annealed in a vacuum and then sealed under vacuum until time 
for assembly. The support system built for the tube consisted of 
a long aluminum beam v;ith six rigid tube supports. This permitted 
the sections to be rotated while being sealed together so that the 
joint was reasonably smooth with no sags or large irregularities. 
After the sections were joined together, the tips of the thermo- 
couple support tubes were broken, the thermocouples inserted and 
sealed with epoxy. The epoxy was then cured for two hours with 
heat lamps. After curing, the tube was placed in its final 
position and sealed to the system via a' section of glass tubing 
containing a metal bellows. 

2 . Pr eparat ion 

Preparation of the tube involved bakeout and outgassing 
of the tube and its elements so that a high vacuum could be ob~ 
tanied followed by filament activation. Bakeout was perfc . ed 
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with heating tapes at 200°C for over 24 hours with the system 

under evacuation. The cold trap was kept filled with liquid 

nitrogen following bakeout . Outgassing was performed by induction 

heating, after bakeout. Since the leads were KOVAR , with a much 

lower melting point than the tantalum cylinders, outgassing of the 

tantalum at elevated temperatures was done in short bursts. In 

the case of the diodes and main discharge electrodes, the filaments 

were activated first and then the anodes were outgassed. Filament 

activation consisted of slowly raising the temperature by varying 

the filament voltage until the temperature, as measured with an 

optical pyrometer, was at a certain value. For the diode filaments, 

this temperature was 1750^C , (subscript B denotes brightness 

B 

temper atui e ) the x ecomiueuded temperature for thcriated filament^. 
Following the manufactur ers instr uct ions , the main cathode fila- 
ments were only heated to 1190^C^. This temperature w^as held for 
five minutes . 

Following outgassing and filament activation, the filaments 
were aged for several hours at their normal operating temperatures. 
Then the main discharge anodes were cleaned by electron bombard- 
ment. This was accomplished by applying 300 volts to the anode, at 
which time the anodes turned bright red under the absorption of 
45 watts of power. 

3. Calibration of the Diodes 

After the tube was prepared, with all elements outgassed, 
the diodes had to be calibrated. Since the mobility-limited cur- 
rent is less than the space-charge limited current, and the 
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transition is a function of gas pressure, it is best first to 
attain space charge limited operation. This was easily determined 
by preparation of a log-log plot of the plate current versus 
plate voltage to a value of plate voltage (30 volts), well in 
excess of that to be used ( 13'i volts). This was found to be a 
straight line, indicating that the filament voltage selected 
provided a good operating point. 

Calibration consisted of varying the pressure in the tube 
and measuring the plate current as the plate voltage was varied. 

The filament voltage was held constant. As the pressure increases, 
there is increased filament cooling. It is not practical to use 
an optical pyrometer to maintain constant filament temperature at 
all pressures. Since the power app3 ied increases as the tempera- 
ture decreases for constant voltage, the constant voltage mode is 
semi -compensatory . 

The filaments were thoriated tungsten and at first operated 
extremely well. After the main discharge was first run, which 
overheated the epoxy seal at the thermocouple, the diode current 
became noisy and started to decrease. It became necessary to 
increase the current to a high value and operate the diode fila- 
ments as pure tungsten. This required a recalibration. 

4 • Data Accumulation 

It was decided to take data at three pressures. The low 
pressure was selected to be around 8 torr because a larger power 
supply was needed to obtain cx good current pulse at pressures much 
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lower than that* The high pressure was l6 torr because of the 
difficulty in getting the main discharge started at higher 
pressures. The third pressure chosen was 12 torr. 

In order to be able to detect the existence of any trends 
(in the data), due to time the order of experiments was chosen 
with a random number table. First the order of the pressures was 
chosen and then the order of the polarity of the main discharge. 
The order of runs is shown in Table II. Due to the difficulty 
of filling to the same pressure, the runs were not completely in- 
dependently or der ed. 



TABLE II 



r\T7^ nnvc 



Polar ity 


A 










L 


+ 1 


0 


-1 


Pr essur e\^ 


L 








8 


1 


2 


1 


3 


12 


3 


3 


2 


1 




2 


1 


3 


2 



A set procedure was followed when taking the final data. 
First all equipment was turned on and allowed to warm up. Then 
the tape recorder was calibrated. The calibration of the oscil- 
loscope was checked. For each pressure, the steps were the same 
except for minor variations due to the different order of the 
polarities within each pressure. The trap in the filling line 
was partially filled with liquid nitrogen. After the tube was 
filled with Argon at the desired pressure, the valve to the filling 
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manifold was secured. The diode and main discharge filaments 
were given plenty of time to warm up. IVlien the main discharge 
was on, the Pupps discharge was adjusted for minimum AC current 
variation in the main discharge. This was observed on the 
oscilloscope using the current probe. The d.c. diode signal was 
measured after adjustment so that the amplified a.c. signal did 
not exceed the clipping level of the slicing circuit. When 
everything was ready, the tape recorder was started. While it 
was running, the temperatures were recorded. Between runs, the 
tape recorder calibration was checked and the output zero adjusted, 
but the input gain was not touched. After the experiment v;as over, 
a 5000 Hz tone burst sine wave, stepped down to 2 millivolts peak 
to peak wa^> feu to the amplifierc in parallel. This wa<= r<:>rorded 
for later use as a calibration aid. After each set of runs at a 
single pressure, the tube was pumped down to a hard vacuum and 
refilled for the next set of runs. 
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V. THEORETICAL ANALYSIS 



This section treats the analysis of data theoretically. It 
develops the expressions needed to utilize the data to obtain 
information from it concerning the processes under investigation. 

In doing so, it takes into consideration the progagation of errors. 

A. ENSEMBLE AVERAGING 

Nine experimental runs were made, each with different conditions. 
An event, or experiment, consists of a current pulse, and the sub- 
sequent sound production. An experimental run consists of many 
events or replications of the experiment. The data recorded for 
each event were the triggering signal and the waveforms at three 
points. After converj^ion Lo digital form, ea^h waveforiu was re- 
presented by a set of samples. Since a run consists of a set of 
events, the data consists of a set of sets of samples which shall 
be referred to as an ensemble of sample sets. 

Consider just one of the signals. The samples in any one of the 
sets were not taken randomly, but have fixed relationships to each 
other which reflects their time dependence. For the set of samples, 
the zero of time is chosen to correspond to the triggering event. 
Thus, the n-th sample represents the value of the signal at time 
t = nT where T is the period of the sampling frequency. 

This sample represents not only the time signal at that time, 
but also something else. It represents the value of the cross 
correlation between the signal at that time and a Dirac delta 
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function signal v/hich coincided with the triggering event. This 
is a subtle point which was used for . advantage in eliminating 
noise . 

The cross correlation function /(x) of two real functions 
g(x) and h(x) is given by: 



/(x) 




CO 

g ( U“X )h ( u ) du 



Let g be the Dirac delta function 6(t) and h be the signal s(t). 

Then * ^ 

/( t) =1 6 (u-t )s (u)du 

-CO 

and therefore, 

/(t) = s(t). 

By ensemble averaging the sets of sarapies , we obtain the 
ensemble average of /(t), which is the cross correlation function 
of the signal with the tr igger ing event. For a large enougli 
ensemble, this should approach the time average of the cross cor- 
relation function. In the time average, contributions from non- 
correlatable portions of the signal, such as those due to noise, 
random or non-random, will vanish. 

Because of their intimate relation, /(t) and s(t) may be used 
interchangeably except that it must be recognized that in each 
sample set, s(t) has a different reference time while /(t) has the 
same reference (for this reason, s(t) will be used for /(t) except 
when the special features of /(t) are needed). 



81 




i 



When /(nT) is ensemble averaged, the mean value and its 
standard deviation are obtained. Since there is no reason to 
assume that the distribution of noise or fluctuations in value 
of /(nT) should depend on n, the standard deviation is considered 
a constant. Then, if /^(nT) represents the value for the i-th 
set of samples, the process of ensemble averaging can be repre- 
sented by 



1 ^ 

/(nT) = - / .(nT) 



( 5 . 1 ) 



N-1 






and 



M 

1 * 

a = — L 



a 



M n 



(5.2) 



where M is the number of points in the sample. It is possible 
to consider this ensemble average /(nT) and the standard deviation 
(7 as the expected value of the signal, s(nT) and its standard 
deviation, cr . 



B. FOURIER TRANSFORMATIONS 



When a harmonic waveform is fed into a linear t ime- invar iant 
device, the output is also harmonic but in general with a dif- 
ferent amplitude and phase [ 48 ]. This is true whether the device 
be a filter, or transducer or some other system with an input and 
an output. Such a device can be completely specified by a fre- 
quency dependent complex quantity Z(v) given by the ratio of the 
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output to the input. If the input v^(t) is given by 

/+\ - A i2nvt 

v^(t) - Ae 



and the output by 

V2(t) 

then 

Z(V) 



„ i2nvT+i0 
Be 

B i0 



If an arbitrary input is applied to the device and it 

is desired to calculate the output v ^{X) , "^1^^^ analyzed into 
its spectrum; each spectral component multiplied by the corre- 
sponding transfer factor to obtain the spectrum of the output, and 
then synthesized from its spectrum. Thus 

V2(V) = Z(V)V^(V) 

and 

V2(t) = ( Z(v)Vj(v)e^^'''^^dv. 

^ —CO 

(Here and elsewhere, where a quantity is viewed in both the 
frequency and time domains, the upper case symbol will be used 
for the frequency domain and its lower case for the time domain.) 
If the process of convolving g(K) and h(x), is represented by 

g*h, defined by ^ 

/(x) = g*h =j g(u-x)h(u)du, 

-'-00 



and the process of Fourier transformation is represented by o , 
the . determination of ^ 2 ^^^ presented by the diagram below. 
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t 



(multiplication) 



v^(t) => V^(V) 



z(t) 

I 



z(v) 

I 



V2(t) 3 V^{V) 



In this experiment, there are three time dependent waveforms, 
i(t), corresponding to the s ound~pr oducing dis- 

charge current, the sound at point 1 and the sound at point 2 
respectively. The object of the experiment was to obtain the 
transfer function of the conversion of current to sound, Z^(v), 
and thef^^^^^G^ function of the spatial region betv^een points 1 
The entire process can be represented by the 
diagram belov/, where the time representation of the transfer 
functions and the convolution process has been omitted for clarity, 



and 2, Z^^{V) 



V 



i(t) ^ i(v) 



X 



So(t) ^ 



Z,(v) 

st(v) 



s^(t) :> 



X 



Zo-,(v) 

i 

S^(v) 



X 



Zi2(V) 






s„(t) S„(V) 
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The application of Fourier transform theory to this problem 
assumes linearity of the processes. That is, if is the 

output at point 2 to an input of i^(t) and if ^22^^^ 
response to i^(t), then the response to i(t) = i^(t)+i 2 (t) is 
s (t)= s (t)+s (t), regardless of the choice of i and i . 

In the transform domain, the process from current to sound at 
point 2 is represented by the set of equations: 



Sq(v) = I(V) 2g(v) 



s^(v) = Sq(v) Zq^(v) 



Since i(t), s (t) and s (t) are known, their respective 

_L ^ 

transforms are known. Therefore these constitute three cquaticnc 



with four unknowns, Sq(v), case 

where the main discharge is on, only \z^^{v)\ is desired. This 
can be solved for easily, 

, , |S2(V)| 



To determine 2 , when the discharge is off, acoustic theory 



provides a model from which a fourth equation can be obtained. For 
sound propagating in a tube, the transfer function between two 
points is given by 



^ 12 ^^^ ^ 



ik(v)d 



12 



where d^^ is the distance between 1 and 2. With this equation and 
the additional assumption that, in the absence of a discharge in 
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the intervening space, h(v) is the same throughout the tube, the 

system reduces to two equations in two unknowns, 

ik(v)x 

S^(V) = I(v) Zg(v)e 

ik(v)x 

S 2 (V) = I(V) Zg(V)e ^ . 

The propagation constant k(v) is complex, 
k(v) = a(V)+ ip (V) . 

Under the exj^er imen tal conditions, boundary losses dominate. In 
this case, P(V) = phase information is neglected, 

meaning that the delay between i(t), S^(t), and S^Ct) is neglected, 
then can be obtained from S^(V), by 




C. ACOUSTIC TRANSFER FUNCTION 

After ensemble averaging, the two diode signals are represented 
by sample sets and standard deviations 

s^(nT), 

n = 1,M 

and s^{nT), 

Since the standard deviations are considered time independent, each 
represents a constant signal in time, the transform of which. S, 
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is a zero frequency component of the same magnitude as 0“. There- 
fore, in the transform domain, the two diode spectra are 



and 

2 2 

where and affect only S(v^ = 0). 

The standard deviation or error in S^(V-) and 

3 . 

to the noise, plus that due to other sources, 
value of a quantity such as is represented 
symbol for the primary quantity, e.g. 

Let Y(v^) represent the natural logorithra 



is that due 
(The error in the 
by A preceeding the 



of Z 



12 



(V.). 



Then 



Y(vp 



= 



Is.{v, )| 



I o / ^ 
I ^1 V 



\ I 
•i i f 



( 5 . 3 ) 



and 




( 5 . 4 ) 



For the discrete case, the earlier expression for is not very 
useful. Another expression which can be written containing is 

Y(vp = - 

This expression can be solved for 6 at each value of v . . However 

o . 1 

a better approach is to use all the information available and 
minimize the error in 6 . If all values of Y(v.) were known to 

the same degree of accuracy, then a simple linear least squares 
fit would provide the optimum estimate for However, from 

the ensemble average and another consideration to be discussed 
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later 5 the error in Y{\)^) is known. A weighted average will 
therefore improve the estimate for The values of are 

best weighted with a weighting function 



W. = 
1 



A 



AY(V. )' 



where A is a normalizing constant given by 



( 5 . 6 ) 



A = S 
i=l 



AY(V. ) 



D. SOURCE TRANSFER FUNCTION 



When there is no discharge between points 1 and 2, it is 
assumed that the transfer function is the same throughout the tube. 
This means that the value of may be used to extrapolate the 



sound spectruiu observed at points 1 and 2 to the crigin. There it 
may be compared with the current spectrum of the sound source 
discharge to obtain the source transfer function. 

The spectral components of S and S^ projected to the origin 
are given by: 



S. (V. ) = S.(V. )g 

jo^ 1 ' J ^ 



P X .V . 

O J 1 



( 5 . 7 ) 



(j=i,2) 

If p^ is used, the implication is that the estimate of p^ is 
independent of the . But since p^ was derived from the estimates 
for S 2 and S^, this is not true. Recall that 





( 5 . 8 ) 
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ther ef or e 



1 






(5.9) 



This can be substituted into Eq. (5.7) to give 



/^l\ >= 2-^1 

= ^(s;j 



and 



S2o('^i) ^2 5 



1 \ 



Since these are equivalent, the estimate for Sq is 



^0 ^2 



( 5 . 10 ) 



The standard deviation for is given by 



AS^ = 



0 0)\X2-X^y 



''2 /iflY ^ / aS 2 \ 2 ) ^2 






.(5.11) 



The transfer function of the source is obtained with a simple 
division, 



I(V.) 



( 5 . 12 ) 



whereby 



2 |''So<'’i>Y As<'’i>V 2 



(5.13) 



E. EFFECT OF THE DISCHARGE 

In order to make a valid comparison betiveen experimental results 
and theoretical predictions, the predictions should be based on 
values which represent the experiment. In this case, to compare 
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the Ingar d-Schulz predictions with experiment, the effects of 
longitudinal and radial thermal gradients must be taken into 
consideration. Since the electron temperature and electron density 
are not known, they must be estimated. Also it should be noted 
that the discharge does not fill the region between the two diodes 
but, using the values from Table 1, occupies only 55-6% of the 
space. When there is no discharge, this does not affect results 
but when the discharge is on, the region is inhoraogenous and the 
attenuation constant obtained is the average value for the entire 
space between the diodes. 

Several assumptions must be made in order to predict the 
experimental results. The first is that there is no transition 



zone at the one 



-'O. iliCiCAiiO i lC\ V C 4 J C. 



predicted value for the effective attenuation constant should be 
the average of the value for no discharge and that in the discharge, 
where the values are v;eighted by the lengths of the two regions. 

Due to the nature of the interface between the regions, reflections 
are considered negligible. Temperature measurements taken between 
the discharge and diode 2 support this assumption. They showed only 
a one-degree higher temperature than measurements far removed from 
the discharge. The two diodes are separated by a distance d^ and 
the main discharge elements by d^. Thus, if and denote 

the attenuation constant in the discharge and outside of the dis- 
charge respectively, the effective attenuations constant is given 

d. 



by 



^0 H f"oi ■ d,. *^02 
0 



0 



where . 
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other assumptions concern the effects of the discharge through 
thermal gradients. During their investigation of the neutral 
temperature in medium pressure discharges, Hayess , et.al. L*49l 
observed a radial temperature gradient under all discharge con- 
ditions. This will affect the propagating sound. Hayess C 5^1 
considered the problem in conjunction with his experiment. He 
concluded that for frequencies below the first radial resonant 
frequency, only the fundamental mode, a plane wave would pro- 
pagate. He further concluded that the temperature difference 
would not affect the attenuation greatly since only its square 
root appears in the attenuation factor. The cut off frequency 
for this experiment is about 9 kHz. The discharge region can 
be charact er ized by two temperatuies for two different purpcccc. 
The longitudinal gradient in the tube affects the pressure in 
the tube. For this purpose, the dischai ge can be char acter ized 
by its mean temperature T^^^. Since the temperature outside of 
the discharge, T^, did not change, the tube pressure with the 
discharge on is given by 




where d^ is the length of the tube* The uncertainty in the 
pressure is given by 
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since 




^9 

2 o 

a. 2 
0 T 

M 



this can be expressed as 




2 o 
d^ T 

0 "m 



2 




2 



Under the assumption that boundary effects are the dominant 

loss mechanism, the temperature at the wall, T^, is important for 

the discharge region. This is because these losses occur within 

two or three mean free paths of the boundary. 

To calculate the predicted attenuation constant with the 

Ingar d-Schulz equations, the experimental parameters must be 

converted to the variables used in the computer program: CO, , 

N , T , and T . Tlie wall temperatuie <^an be used for T . The 
e e n ' n 

neutral density can be calculated from p’ with T^ , using the 
expr ess ion 

1 

N = 9.6642 X 10 p/T 
n ^ 

-3 . . . 

where p is in torr and N in cm . The relative uncertainty in 
^ n 

the density is given by 




Since they were not measured, estimates must be used for the 
electron temperature and the electron density. With the similarity 
parameter Rp as an entry, T^ may be obtained from von Engle [4o] 
with an estimated uncertainty of 10%. Estimating the electron 
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density requires a more complicated process. Venzke C51] has 
plots which give E/p as a function of i/R with Rp as a curve 
parameter. The value of E/p thus obtained can be used to obtain 
the electron drift velocity from page 5^5 of MdDaniel [39]. With 
the assumption of a uniform current density across the tube, the 
electron density is given by 



N 

e 



2 nr V q 
o D^e 



Of course it is known that the current density is not uniform 
across the tube but neither are other factors such as neutral 
density. The* entire acoustic problem is treated in a one-dimensional 
sense. Because of the approximations involved, the electron den- 
sity is assumed to be uncertain by 50^. 

The relation of the attenuation constant, to p, the 

attenuation co-efficient, can be viewed in two different ways. 

If p = p^ , the two ways are equivalent. One way is to write 

PQ = P(W = 1). 

The other way is to write 

Po - P(^)/ 

The latter form is preferred since a departure from the assumption 
of boundary loss domination will be readily apparent from a plot 
of p^ as a function of U). It should be noted that Figure 2 , which 



is a plot of 



Pi-Pui 

^UI 



, v/ould not be changed by switching p for p 



Since p^ should have no dependence on CO, the value of co should not 
be important. To be certain, the value corresponding to the highest 
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and lowest e>:per imental values of 60 should be used and the 
dependence checked . 

Since all of the input quantities are estimates, with an 
uncertainty which can be treated as a standard deviation, the 
estimate of obtained will be uncertain. The uncertainty in 
is given by 




The partial derivatives can be obtained with the computer program 
by the expression 

be e (x.=X.+AX.) - B (x.=X.) 

bx. AX. 

1 1 

If Ap^ is particularly sensitive to any of the variables, the 
process of evaluation should be reviewed to see if the uncertainty 
can be reduced. 
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VI, TREATMENT OF D\TA 



The choice of data handling procedure was strongly influenced 
by problems associated with the experiment. The strong reflection 
from the tube ends was one of the most important factors in all 
considerations. The need for gated signals rather than continuous 
waves, in order to study the effect of the discharge, meant that 
the data could not be taken with meters. The multi-component 
spectrum associated with the bias for the gated sine waves meant 
that the signal had to be Fourier analyzed. 

Due to the strong reflections, the signals could not be 
analyzed in the continuous mode with a Correlation Function 
Computer such as the PAR lOOA . When the gated mode was used, the 
time base was limited to the period from the triggering signal to 
the start of the reflected signcxl. Since this time base was 
sampled only 100 times, the resolution was very poor. Preliminary 
analysis of the signals indicated that the source transfer function 
decreased at low frequencies. Since this could help to check 
theories of the production mechanism, the low frequency components 
were of particular interest and resolution was needed. 

A Boxcar integrator was also used to recover the signal on a 
scale compatible with transcription into numbers. However like 
the Correlator, this utilized only one signal waveform at a time 
and it wasted most of that. In addition, the data had to be con- 
verted into numbers, and card punched for computer processing. 
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Even using digital voltmeters to obtain the data, obtaining more 
than one hundred twelve points from a waveform was very difficult. 

For these reasons, the process of recording on tape and 
direct conversion to digital form was finally selected. It offered 
the greatest resolution, processing possibilities and error re- 
cover y capability . 

A. CATA REDUCTION PROCEDURE 

1 . Analog to Digital Conversion 

After completing the data recording, the data tapes were 
reviewed with an oscilloscope. One minute segments were cut from 
each of the original tapes and spliced together. They v/ere 
separated by long sections of white leader tape for easy identifi- 
cation of the segments, 

A hybrid computer composed of a COMCOR 5000 analog computer 
interfaced with an SDS 9300 digital computer was used for the con- 
version of the analog data to digital form. The analog side of the 
hybrid computer was set up to receive the analog signal and provide 
the appropriate signals to the digital portion for digitizing the 
information. The gain of the data channels was 50. A channel 
corresponds to the waveform recorded on it. Logic boards were 
used to start sampling at a frequency of 10 kHz w^hen the trig- 
ger signal reached a variable trigger level. Two switches were 
used to start and stop the conversion process. The number of 
samples taken following each trigger signal v;as controlled by the 
digital program. After converting 1200 samples of each of the 
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three data channels to digital form, the information was written 
on a magnetic tape. The program used for the conversion is a 
standard program available to all users of the hybrid computer. 

The parameters for the program are the number of channels to be 
sampled and the number of samples per channel to be made for each 
record. The logic involved in the sampling process is different 
for most conversions and this required a slight amount of patching 
on the logic patch boards. 

2 . Seven to Nine Trac k Conver s ions 

The information written on the 7 track magnetic tape by 
the hybrid computer is in the form of 2^ bit binary words. To be 
processed by the IBM 3^0, these must be converted to 32 bit binary 
words. The only programming involved in lijis is a prograia to read 
the 7 track tape, call a conversion routine and write the information 
on another tape, usually 9 track for increased efficiency, for later 
processing. For this job, this step was also used to separate the 
samples by channel. 

3 . Ensemble Averaging with Data Rejection 

The diode signals contained noise. The noise was random and 
consisted chiefly of positive-going spikes. A low pass filter could 
not eliminate them without distorting a significant portion of the 
acoustic waveform. In addition there were a few bad trigger signals 
recorded. Thus it was known from the beginning that some of the 
data records were bad and should not be included in the processing 
but it was not known which they were. Since the entire data set 
consisted of ten files, each of which contained from 80 to . 0 
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records of 36OO words, screening became a major task. The 
following procedure v;as finally found to be most efficient and 
effective . 

Each file corresponded to an experimental run and all 
records in that file represented the same waveforms plus noise. 

The records in a file were all averaged together and the standard 
deviation for each channel was obtained. Following this, all 
records were compared with the average. Those with less than 75 % 
of their points within 2.8 standard deviations of the average were 
dropped from further processing. Using only the acceptable records, 
a new ensemble average and standard deviation was obtained for 
each channel. These were used in a final ensemble averaging pro- 
cess in which only those Doinxs within 2 standaru dev/iatious of 
the average were used. This provided an ensemble average of the 
waveform and the standard deviation of the points involved. 

4 • Delay, Overlay and Zero 

At this point, the data consisted of 1200 samples for each 
of the three waveforms, representing their average values, and the 
standard deviation of the waveforms. For the sampling process 
1200 samples were taken so that all of the desired waveform of the 
signal at point 2 was sampled. This meant an undesired signal was 
included in the samples of the signal at 1, namely the reflected 
signal. In order to make a valid comparison of spectral components, 
the time samples of the waveform which is transformed must contain 
the same number of samples. Reducing the number of samples of to 
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those necessary to avoid the inflections in s^ means the starting 
point for the transformation must be delayed. Since phase infor- 
mation Vv^as not going to be considered in the analysis, and time 
delay merely changes the phase of the components, delaying the 
signals would not affect the results. Therefore the signals s^ 
and s^ were delayed by choosing a starting point for each such 
that 600 samples included the waveform up to the point where the 
reflection started. 

The diodes interacted electrostatically with discharge. 

This caused a small coherent square v;ave in the diode output 
which coincided with the time of the current pulse. Its’ inclusion 
in the waveform to be transformed would give erroneous spectral 
components for the acoustic signal. Delaying was sufficient to 
eliminate this unwanted signal from s^. In order to eliminate it 
from s^, without introducing unwanted spectral components, the 
values of the samples in the region containing this unwanted 
signal were replaced with values taken at random from the adjacent 
region. Since this adjacent region contained only samples of the 
waveform before arrival of the acoustic signal, they only repre- 
sented background noise. 

Since the initial portions of the signals S^ and S^ 
represent no signal, i.e. the period before arrival of the acoustic 
signal from the sound source, its average value should be zero. 

The signals were offset to zero by first determining the average 
value of the no signal region and then subtracting it from all 
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values. For the current waveform, this region followed the 
current pulse. The three operations of delaying the signals, 
overlaying the electrostatic signal and establishing the zero 
level were performed in one program, 

B. OPERATIONS ON THE TRANSFORM 

Two sequential programs were used to transform the vjaveforms 
and obtain acoustic informations from them. The first program 
took the input waveforms, corrected them for gain and diode 
sensitivity and determined their spectra. The second used the 
spectra to determine the transfer functions of the discharge 
region and the source. 

Input data for the transform program were the three ensemble-' 
averaged waveforms, i(nT), s^(nT), S 2 (nT) and their standard 
deviations, /^i, AS^ j ^iid the diode sensitivities. 

Two uncalibrated amplifiers were used in amplifying and 
recording the diode signals. These were calibrated as part of the 
processing. While recording data, both amplifiers were fed a 
common signal consisting of 2 millivolt peak to peak, 5000 Hz, 
tone burst. This calibration signal was processed in the same 
manner as the other signals. In the first part of the transform 
program, the spectrum of this signal, in the diode channels, was 
Fourier analyzed. Then the spectral component closest to 5000 Hz 
(5040) was used to calibrate the gain for each diode channels so 
that this component was equal to one. 
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The diode sensitivities, supplied as inputs, were calculated 
earlier from experimental data and expressed in terms of milli- 
volts/torr. The sensitivity of the current detector was 
5 X 10 volt s/mill iamp . The sensitivities and corrected gain 

factors were used to convert the values into the correct magni- 
tudes of the physical quantities involved, milliaraps for the 
current and torr for the acoustic signals. 

The transforms of the corrected waveforms were then calcu- 
lated using a computer library subroutine '^RHARM^\ This program 
requires that the number of points be an integer power of two. 

There were 600 points available and the nearest power of two was 
512, leaving 88 excess points. This means that there is some 
u n c er t a i i j t y as to w I'l er e to start the an alysi-S « Per this reason , 
the waveforms were analyzed five times, each time advancing the 
starting position by 10 places. The average value and variance 
of each component was calculated from the five spectra. 

The output of the transform program was the average spectrum 
and the variance of its components. Only the first ‘twenty five 
of these were passed to the next program. The magnitudes of the 
remaining components were less than 5% of the fundamental and were 
characterized by a great deal of scatter. 

The next program used this data, and calculated the acoustic and 
source transfer functions. The equations of Chapter V were used 
for the calculations. A computer library subroutine ”LSQPL2^^ was 
used to obtain the least squares determination of For this 



101 



calculation, the data was weighted inversely proportionally to 
the Square of its uncertainty, with weights normalized such that 
their sum was equal to the number of points. 

For the data corresponding to no discharge, Eq. (5*10) was 
used to obtain an estimate for the sound spectra at the source. 
This was divided by the sound source current spectra to obtain 
the transfer function of the sound source. 
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VII. PRESENTATION AND DISCUSSION OF RESULTS 



A. AVERAGE WAVEFORMS 

The averaging and screening method worked very well although 
the program could be rewritten to perform the job more efficiently. 
The three waveforms for the discharge-off case at 7*65 torr are 
shown as Figures 10, 11 and 12 just before being transformed. For 
contrast, the first 2/3 of the average \;aveform for s^ for the 
same run is shown as Figure 13 with the same scale. The second 
pulse is the reflected signal from the end close to the sound 
source. Thi.s shows the effect of the delaying, overlay and zeroing 
operations. Unless otherwise specified, the remainder of the 
figures displayed in this chapter were derived from these wave- 
forms . 

B. ACOUSTIC TRANSFER FUNCTION 

The points of the acoustic transfer function showed a rea- 
sonable amount of scatter about the predicted form. The greatest 
degree of scatter was found to be associated with those components 
whose contributions to the signal were small. Figure l4 shows 
the first twenty five spectral components of the sound at point 2. 
Figure 15 shows the transfer function in a linear form. Figure l6 
is the same information expressed as attenuation of sound in 
decibels per meter as a function of frequency. The curve repre- 
sents the least squares fit to the data under the assumption that 
boundary losses dominate. The goodness of fit found in all cases 
supports the validity of this assumption. 
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Figure 10. Sound Source Discharge Current. P 
Main Discharge off. Time Scale: 2 
Vertical scale is arbitrary. 




Figure 1]. Sound Signal at Point 1. (delayed, 
and zeroed), P = 7.65 torr. Main 
off. Time scale: 2.048 msec/in. 
scale is arbitrary. 
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Figure 12. Sound Signal at Point 2, (delayed, ovei layed 
and zeroed). P = 7.65 torr . Main discharge 
off. Time scale: 2.048 msec/in. Vertical 
scale is arbitrary. 




Figure 13. Sound Signal at Point 1. ( Uncoi r ected ) . 

P = 7.65 torr. Main discharge off. Time 
scale = 2.048 msec/in. Vertical scale is 
the same as Figure 11. 
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P'igur G l4 . Si'-ectral Decomposition of Sound Signal at I c t 
TJ'.is is the -spectral c’eccnpcs j t ion C'f , ■ ■ ■ 
form shown in Figure Id. Dashed line -- - 
simple connection of points, 
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Figure 15 . S^{v) /S^{v) in Dependence on V ^ for a Pressure 
of 7.65 torr . Main discharge off. Error bars 
represent one standard deviation. Curve 
represents the weighted least squares fit to 
the data. 
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Figure l6. Attenuation in Dependence on Frequency for a 
pressure of 7.65 torr . Main discharge off. 
Curve represents the least squares fit to the 
data. This figure corresponds to Figure 15. 
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C. SOURCE TRANSFER FUNCTION 



The source transfer function was found to have the low fre- 
quency drop off indicated by the preliminary observations. At 
:ffe<^uencies above a pressure-dependent transition frequency, the 
frequency characteristic did follow the predicted form. The 
magnitude was also found to agree with the predicted values. The 
transfer functions are displayed as Figures 17, l8 and 19 for the 
filling pressures 7-65 torr , 12.7 torr and l6 torr respectively. 

The structure is believed to be due to scatter in the data since 
the points with the greatest departure from the smooth curve 
correspond to those spectral components with the largest relative 
err or . 

1. Agreement with Pre di ction 

Two aspects of the results can be used to test the theory 
of sound production from a discharge. One is the frequency de- 
pendence and the other is the magnitude, 
a. Frequency Dependence 

The low-frequency drop off is contrary to the fre- 
quency dependence predicted by Schulz's results. This effect is 
believed to be due to neutral particle diffusion under the in- 
fluence of thermal gradients. It is a low-frequency, lov;-pr essur e 
effect which should appear only when thermal energy sources provide 
the acoustic energy. This will be discussed later. 

According to Schulz's prediction, the transfer function 
should vary as 

CO 
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Figure 17- Source Transfer Function in Dependence on 
Frequency for p = 7*65 torr. Main dischar 
off. Vertical scale in units of (newton/ri , 
per ampere. Solid curve represents predicted 
form fitted at one point.. 
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Figure l8. Source Transfer Function in Dependence on 

Frequency for p = 12.7 torr . Main discharg 
off. Vertical scale in units of (nevvtons/m ; 
per ampere. Solid curve represents predicted 
form fitted at one point. 
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Figure 19. Source Transfer Function in Dependence on 

Frequency for p = l6 torr , Main discharge off. 
Vertical scale in units of (newton/m^) per 
ampere. Solid curve represents predicted form 
fitted at one point. 
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This prediction was tested by fitting a curve given by 



1 . , 0 &, 
y = TT- sin ( — ) 



(7.1) 



through the transfer function data points * The coefficient 
was evaluated by causing the curve to pass through the point at 
V = 2450 Hz . 



Schulz’s derivation was based on planar electrodes 
separated a distance ’’a^\ The experiment used cylindrical elec- 
trodes of length 1^ cm separated by a distance of 3 cm. Therefore 
the effective length of the discharge could lie anywhere between 
3 and 6 cm. The value of ^^a’’ was adjusted within these limits to 
obtain a curve which fitted the data reasonably well. For 7-65 
tort, the value of a was 3-5 cm: for I 6 torr , the value was 
4.6 cm. Since the pressure affects the discharge, there was a 
possibility that the value of ^^a" for 12.7 torr could be inter- 
polated from the other two values. The value obtained by inter- 
polation was 3*93 • This gave a good fit to the data. The smooth 
curves shown in the figures are the fits of equation (7.1) to the 
data . 



b. Magnitude 



Using Schulz’s equation, the source transfer function 



is predicted as 



Zg = (y-i)H 



ICO 



The energy transfer rate from the electrons H is the 
same as used earlier in conjunction with wave amplification. It 
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may be expressed in terms of the characteristic time for electron 



heating as 



H = 



T(y-l) • 



The calculated value of c /T v;as obtained with the 

computer program used for evaluating the Ingar d-Schulz equations. 

This was accomplished for the transfer function at v = 2540 Hz 

and p = l6 torr . The discharge current spectral component at 

this frequency was 6.4 ma which was estimated to correspond to 

10 -3 

an electron density of 2 x 10 cm . In determining this density, 
the cross sectional area of the electrodes was used since the dis- 
charge did not extend, in the radial direction, significantly be- 
yond them. The procedure used was the same as outlined in Sec. V.E. 
The predicted value of obtained was 

^ S , 2 

r- nt/m ^ ^ .. - _ . 

-LW U-/X JL J 



_ , ^ r' nt/m 

(M -V 'J ■■■ ■ - ■! II ^ H II 

"S * ampere 



= 43.5 • (V = 2540 Hz, p = 



The corresponding value, experimentally obtained, was 

Z = 33 2 
S ampere ‘ 



The agreement of the two values is good. Considering 
the number of approximations made, order of magnitude agreement 
would have been considered satisfactory. No error analysis was 
conducted on this computation for that reason. 

Except for the departure at low frequencies, the re- 
sults are consistent with the theoretical predictions based on the 
assumption of energy transfer by the electrons in the positive 
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column . There are aspects of the problem which were not examined 
by this experiment but which should be studied to gain complete 
understanding and to understand other e>rper imental results. A 
more complete investigation would include phase as well as magni- 
tude information in the transfer function. In addition, the 
source of the sound could be localized and momentum transfer 
aspects investigated. Although this experiment indicates that 
electron heating of the neutrals in the positive column provides 
a mechanism for sound production in glow discharges, it has not 
discounted the roll of ions in sound production under other 
circumstances . 

2 • Low Frequency Departure 

During preliminary work in conjunction with this c:q^cri- 
ment , a large number of acoustic waveforms were observed. One 
characteristic of the sound signal produced by a pulse v;as noticed, 
which is reflected in the drop-off of the transfer function at 
low frequencies. Figure 20 shows part of the waveform received at 
point 1 during one of the runs. EXie to attenuation of the high 
frequency components, the features are not as sharp as some which 
were observed at closer positions during earlier work. The main 
features are visible however. 

The first and most significant feature is that the vertical 
distance from 1 to 2 is approximately equal to the distance from 
3 to 4) with the latter always slightly greater. The next feature 
is that the slope and curvature of the curve from 2 to 3 is approxi- 
mately the negative of that for the curve from 4 to 5* Th t vo 
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Figure 20. 



A Typical Sound Pulse Produced by a 
Rectangular Current Pulse. 
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curves suggest an exponential process with a long time constant. 

A last point is that point 5> vvhere the slope is approximately 
zero, is lower than point 1, the level before the signal. The 
polarity of the figure is such that point 1 corresponds to the 
ambient pressure with upward deflections indicating an increase 
in pressure in the discharge region. 

The following description of events occurring in the dis- 
charge region during the creation of the sound pulse is postulated. 
When the discharge current increases, the energy transferred to 
the neutral gas increases, thereby increasing the temperature of 
the neutral particles. This temperature increase is so rapid that 
the netural particle density remains essentially constant. As a 
consequence, the pressure changes almost directly w^th tem- 

perature. The neutral particles move away from the discharge, 
radially as well as axially, under the temperature and pressure 
gradients. As the concentration decreases, the pressure decreases 
toward the ambient value. The pressure change, Ap^, which would 
take place because of the diffusion is 

Ap^ = KANT^ 

where AN is the change in density. 

When the heating processes are removed, the neutral gas 
cools rapidly to the ambient temperature- Since the cooling 
process is much more rapid than the particle diffusion, the pres- 
sure at ambient temperature would be the equilibrium pressure 
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minus the pressure difference caused by the density change. This 
pressure difference, given by 

Ap^ = KANT^. 

Since Ap^ ^ neutral particles would then 

diffuse back in, to return the region to the ambient density. 

The diffusion process is slow, with a long time constant. 
Under the action of a high frequency signal, there would not be 
time for a significant amount to occur between cycles. Consequently, 
the acoustic pressure would not reflect any density changes. At 
the other extreme, for very low frequency temperature variations, 
the pressure changes would be negligible but the density changes 
would be large. Therefore it can be predicted that the effect 
will become appreciable at low frequencies. 

Notice that in the postulated process, no assumptions were 
made concerning the heating process except that it is spatially 
uniform in one region and that particles are free to move in and 
out of the region. Thus the assumption of one dimensionality for 
such a process occurring in a tube could be a valid assumption. 
Therefore Schulz’s assumptions will be considered valid and the 
failure of the prediction at low frequency attributed to the 
wave equation. The wave equation fails because it does not take 
diffusion into consideration. 

If diffusion is the important factor in the discrepancy at 
low frequencies, then the effects will become important for fre- 
quencies whose period is of the order of or longer than a character- 
istic diffusion time T . This time is the time for partic-.e 
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diffusion for a mean displacement of some characteristic dimension 
of the discharge region. Therefore, the characteristic frequency 
is given by McDaniel C 39 ] as 



V = i 

c - T - 

where C is the characteristic dimension. The effective diffusion 
coefficient D is given by Present C36] as D = where v is the 

mean velocity and i is the mean free path. 

Since the mean free path varies inversely with pressure, 

this predicts that the frequency at which the effect occurs will 

vary inversely with pressure. An examination of Figures 15 , I6 

and 17 shows that this is the case: at I6 torr , ^ 500 Hz; at 

12,7 torr, V ^ 750 and at 7.^5 torr, v ^ IGOO \lz . the product 
c c 

pV^ being roughly constant as it would for the predicted dependence 
on density. However, substitution of the discharge length ”a’* for 
the characteristic dimension gives a frequency which is three ordeis 
of magnitude too low. Thus, although diffusion effects may be 
responsible for the low frequency effects, the simple argument 
used to determine the critical frequency did not predict the 
correct result. 

This diffusion effect is a new problem for acoustics. It 
is a low frequency, low pressure effect and will probably only 
occur in conjuction with a thermal acoustic source. The excitation 
of sound by a heat or energy source is self -rectify ing since the 
energy transfer into the system cannot be balanced by the removal 
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of heat or energy from the system by the same process. This 
means that the energy source term H in equation (2.1) can not be 
represented by the form 

iCOt 

H(r,t) = H(r)e^“^ 

since such a representation implies that energy can be removed by 
the source with the same efficiency that it transfers energy into 
the source. 

To predict the correct low frequency characteristics of 
sound produced by a discharge plasma, the derivation of the wave 
equation must take diffusion and heating effects into consideration. 
The results of the derivation should be a differential equation 
capable of predicting the neutral density and temperature when 
solved for the zero frequency case. Such a derivation is con- 
sidered beyond the scope of this investigation. 

D. EFFECT OF THE DISCHARGE 

1 . Predicted Attenuation 

When the measured values of axial temperature v/ere com- 
pared with the predictions of Venzke et.al. CSllj the agreement 
was within 7 %- This verification of the predictions of that paper 
was considered sufficient justification to use its predictions for 
T^^, the wall temperature and T^^, the mean temperature with the 
same relative error. These values were substituted into the ex- 
pressions, of Section V.E. to obtain the neutral density in the two 
regions for each run. The steps outlined in that section v;ere 
followed to obtain estimates for N and T . With the density, 
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wall temperature, electron temperature and electron density as 
inputs, the attenuation constant was calculated using the Ingard- 
Schulz equations for CO = 31j 400 and CO = 600 . 

The results obtained were that the zero ionization value 
of the attenuation constant was essentially the same as the value 
predicted for a discharge under the experimental conditions. The 
value of calculated in the manner mentioned earlier was larger 
at the higher frequency, which indicated that bulk losses were not 
completely negligible. The largest difference between the values 
of for the two frequencies occurred at 7-65 torr where the 
difference was slightly less than 3% ’the un-ionzed value. The 

greatest ionization effect predicted was a , 3 % reduction in 
attenuation for the tv:o discharge runs made with a filling 
sure of 7.65 torr. In other words, the Ingar d-Schulz equations 
predicted that no plasma effects would be discernable under the 
experimental conditions . 

For the plasma associated with a filling pressure of 

7.65 torr, the electron temperature T^ was estimated to be 

1.4 X lO^K and the electron density was 5 x 10^^ cm With these 

17 -3 

values and the values for the neutrals, N = 2.4 x 10 cm and 

n 

T^ = 323 K, the values of some of the plasma parameters were cal- 

9 

culated. The characteristic collision frequencies are: *60 = 2.9 xlO , 

^ en 

7 *"1 

and CO. = (jO = 5*87 x 10 sec . The electron and ion plasma 
in nn ^ 

frequencies 4.00 x 10^^ and 1.5 x 10^ cm The mean 
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cm. While these 



.Zi 

free path of a neutral particle is 7.0 >: 10 
values are different for the other discharges, the values are 
representative of the orders of magnitude of the various para- 
meter s . 



2 . Comparison of Predicted and Experimental Values 

The cxssumption of boundary loss domination was made 
earlier. Since the failure of this assumption would occur at 
higher frequencies, where the data were weighted the least due 
to their larger relat iveecr or , the fit to the data should reflect 
the low frequency characteristics. This was tested at the value 
for predicted for O) = 600. The predicted values are shown in 
Table III under the column marked The experimentally deter- 
mined value of P and its standard shown nn tho 

* o 

columns marked p and A6 . 

The ratio of the difference, between the predicted and 
experimental value of p , to the standard deviation of p is 

° Po-p; 

shown in the column headed r . The magnitude of these 

AP 

*^0 

numbers are all less than unity, indicating that the differences 
between the predicted and the experimental value are all within 
one standard deviation of the result. The greatest difference 
between the predicted and experimental values is 13 %. 

However the signs of the difference are correlated. The 
second column, headed "polarity", indicates the polarity of the 
discharge with respect to sound propagation. A plus sign means 
sound propagation in the direction of electron drift motion. A 
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COMPARISON OF PREDICTED AND EXPERIMENTAL VALUES OF ATTENUATION CONSTANT 
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Uncertainty in predicted value of due to uncertainty in values of wall temperature 
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zero is used to indicate the discharge-off state and a minus sign 
indicates sound propagation counter to electron drift. The sign 
of the differences is positive for the no-discharge cases. For 
the discharge cases, it is negative. The negative sign indicates 
that the predicted attenuation is greater than the observed 
attenuation . 

3 . ■ Discussion of the Results of the Comparison 

The results show that for the ionized cases, the predicted 
value of the attenuation constant is greater than the observed 
value. Although the two values lie within one standard deviation 
of the experimental value, the consistency of the result indicates 
that it is due to either a systematic error or to a plasma effect. 
The temperature dependence of the attenua+ion constant at ronstant 
pressure is such that if the vjall temperature used for the com- 
putation were too low, the predicted values of would be too 
high. The three parameters which influence the predicted values 

to any degree are CO, N and T . The error must be due to either 

n n 

N or T since although the value of CO is the same for the ionized 
n n 

and un-ionized cases, the results of the comparison show that the 
sense of the errors involved is different. 

In calculating the predicted values, the propagation of 
error v.^s followed. Although the mean temperature in the dis- 
charge was assumed to be uncertain by 10%, it did not affect the 
error of the corrected pressure significantly. That error was 
affected primarily by the uncertainty in the value of the filling 
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pressure, which was estimated to be 5%. If the error in filling 
pressure was systematic, it would have affected all cases in the 
same manner. Therefore neither an error in the mean temperature 
nor an error in the filling pressure is believed to be a source 
of systematic error. 

The only remaining source of systematic error would be the 
value of wall temperature used. This affects the results through 
the density, calculated using T^ and p^and directly. The error 
in the wall temperature was conservatively over-estimated to be 
10^, which means a 25% error with respect to the increase above 
ambient temperature but only a 3% error with respect to the ab- 
.solute temperature. Representing the error in due to the 
error xii T b y j-o given by 






I 5n 

\ n 



r 

5T 



AT 



ii) 



Since 



for constant pressure. 



N, 






KT 



5N _ -N 
5t T ’ 



which gives 

*^T ^ T bN^ bT y 00 

Using this expression, was calculated for = 10. The 

results are displayed in Table III. It is apparent that these 
values are not large enough, to account for the discrepancy between 
results except for runs 7 and 9. It should be noted that if the 
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value of the filling pressure used for predicting for runs 7, 

8 and 9, were high, then the correct value for the prediction would 
be higher for all three cases. This would decrease the dis- 
crepancy for the no discharge case and cause the results for the 
ionized cases to be more like the others. The conclusions to 
be drawn from this discussion are that while a systematic error in 
the value of the wall temperature used for predicting the results 
is possible, it is not likely to be the cause of the difference 
between the experimental and predicted results. 

The two runs where this error should be the smallest are 
runs 2 and since the temperature increase is less at lower 
pressures. The difference for these two runs is not significantly 
lower. If the results for filling pressures 7-65 and l6 torr alcnc 
are compared, then it appears that the polarity does have an effect. 
However, there is insufficient replication to support such a con- 
clus ion . 

The conclusion supported by the results is that for these 
experimental conditions, amplification of the sound" may have 
occurred and a polarity effect may be present. The observed effects 
are within the experimental error. 

E. EVALUATION OF THE YATSUI EXPERIMENT 

Yatsui, et.al. [ 8 , 9 ] conducted their experiment in Argon at 
low pressures and high frequencies. For a comparison of their 
experimental results with the Ingar d-Schulz prediction, the 
Ingar d“Schulz equations have now been evaluated with parameter 
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values corresponding to two points in the Yatsui experiment. For 
V = 10 ^ Hz and 1 = 5 niA, they reported maximum amplification 

r- -2 

occurred at 1.5 x 10 torr and damping occurred for pressures 
greater than .31 torr. For the condition of maximum amplification, 
the ratio of the received signals with and without a discharge was 
3.3* The predicted value for this set of experimental parameters 
will be discussed first. 

Since the signal traveled only a fraction C of its total 
travel in the discharge, the ratio of the two signals is 



v' -(ep +(l-e)p)d -£(p -p)d 

S-p _ ^ 



V, 



.pd 



S e 

where d is the total distance between source and receiver. There- 
fore, the difference, between the two attenuation coefficients 

can be determined from their data by 



B -B = 071 

Pp P ed 



The denominator, ed, is the distance traveled in the discharge, 
38 cm. Therefore the experimentally determined value is 

.-2 -1 




-1 

38 



p -p = 3.3 = - 3.12 X 10 cm 



Using the procedures outlined in Section V.E. , the electron 

4 

temperature was estimated to be 5.5 x 10 K and the electron 

10 -3 

density was estimated to be 5.1 x 10 cm . With these values, 
and U) = 6.28 X 10^, ~ 5 x 10^^ cm and T^ = 300K , the quantity 
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p -p was computed using the computer program. This calculation 
gave Pp“P = - 5*3 x 10 , which is lower than the experimental 

value by two orders of magnitude. Using liberal estimates of the 

errors in the input variables, the uncertainty of this result is 

-3 

A(pp-p) = 1.52 X 10 . This is not enough to give order of 

magnitude agreement. 

The Ingar d-Schulz predictions were also evaluated for the 



discharge conditions corresponding to the critical pressure, 

17 10 

.31 torr, using: N = 1.1 x 10 cm , N = 5-5 x 10 , 

n . e 

T^ = 1.6 X lO^K and T^ = 3 OOK . The value predicted was found to 

t>e r. A rff “I 

P -ft = - 9.54 X 10 cm 
P 



Since this difference would hardly be detectable, it is in 
agreement with the Yatsui results. Further agreement with the 
results is the prediction that, for pr essures ‘ above this, damping 
occur s . 

Thus, although the Ingar d-Schulz equations do not predict a 
directional effect, its predictions are, to a certain degree, in 
agreement with experiment. 
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VIII. SUMMARY AND CONCLUSIONS 



This study has examined the creation and modification of 
sound by a discharge plasma both theoretically and experimentally. 
Following a review of earlier experiments on the creation of sound, 
a mechanism for sound production by ions was proposed and shown 
to be feasible. Then a numerical analysis of the equations of 
Ingard and Schulz was performed. This analysis showed that net 
amplification would not occur. It also showed that the approxi- 
mate form of the dispersion relation was not useful because the 
criterion for its applicability was invalid. 

Incidental to the experimental study, the theory of the 
mobility-limited thermionic diode used as a microphone was de- 
veloped beyond its previous, elementary form. This analysis and 
the subsequent successful usage as an acoustically transparent 
microphone for two point measurements has increased its usefulness 
as an experimental tool. 

The diodes were used in an experimental investigation of the 
modification and creation of sound at medium pressures. Fourier 
transforms and computer data processing techniques were used in 
the analysis of the data. The investigation into the modification 
effect of the discharge was in a different region of current, 
pressure and frequency from earlier experiments. The results in- 
dicate that an amplification mechanism may be operative and that 
there may be directional effects. Since the Ingar d-Schulz 
equations, evaluated for the experimental conditions, predicted 
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no effect would be detectable, these results do not support the 
equations. The Ingar d-Schulz equations were also evaluated for 
the conditions of the Yatsui experiment. The predicted values 
were not found to be in agreement with the experimental results 
but the trend with pressure was in agreement. The conclusion 
drawn was that the theory needed more development. 

The major contribution of this study is believed to be the 
information gained from what was originally the secondary experi- 
ment, the conversion of discharge current fluctuations into sound. 
This is the first known investigation in which the current magni- 
tude was compared to the sound it produced. The results were in 
reasonable agreement with the Ingar d formulation of the problem. 

A low frequency effect was observed which was interpreted as the 
first observation of the effects of particle diffusion in acoustics. 

The general agreement of Ingards predictions in both areas 
indicates that electron heating of the neutrals is the dominant 
source of acoustic energy. However several experiments have 
raised questions not answered with the present theory. One of 
these is the directional effect on discharge modification. Two 
others, associated with sound production, are the effects of high 
fields and the localization of the sound source near the cathode. 
These stand out as subjects of further experimental and theoretical 
development. The interpretation of the low frequency behavior of 
the sound source as a diffusion effect also warrants further work, 
both experimental and theoretical. In addition, the study of the 
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Ingar d- Schulz equations indicates that the increased attenuation 
to be expected under certain conditions could be a fruitful area 
for investigation . 
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APPENDIX A 



THE MOBILITY-LIMITED THERMIONIC DIODE AS A MICROPHONE 
Standard microphones have limitations which hinder their 
usage when studying sound at very low pressures. This is parti- 
cularly true in conjunction with plasma-acoustic work. There are 
several reasons for this. A very important one is contamination 
from gassy components. Another consideration is that the micro- 
phones were designed for operation at atmospheric pressure. 
Operation at reduced pressures may require circuit modification 
since, for condenser microphones, breakdown can occur. Chalupnik 
et.al. C 52 J reported changes in sensitivity and frequency response 
with pressure in two commonly used condenser microphones (Western 
Electric 6ZfOAA, Bruel and Kjaer type 4lH). The lack of a sensi- 
tive sound detector without the problems outlined above has 
slowed plasma-acoustic studies. Recently Dayton et.al. [ 53l de- 
ported the use of a mobility-limited diode as a microphone without 
the drawbacks of conventional microphones at low pressures. How- 
ever heretofore the theory has not been developed beyond its most 
elementary form. 

1. Theory 

If the cathode is sufficently hot, so that the operation 
is not emission limited, the current flow in a thermionic diode, 
operating in a high vacuum, is limited by the effects of space 
charge. The current density is given by the La.ngmuir-Childs law 
which says that the density is proportional to the cathode-plate 
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potential to the three-halves power. The constant of pro- 
portionality is dependent on the geomentry of the current-flow 
situation. The three halves power dependence assumes an equi- 
potential filament; a condition which is approximately true pro- 
vided the voltage drop across the filament is small compared to 
the cathode-anode potential. If this condition is not met, then 
the current density varies with the plate potential to a higher 
power, up to 5/2 [54J, the actual value of which is determined by 
the definition of the cathode to anode potential, 
a. Mobility-limited Operation 

In the presence of a background gas of neutral parti- 
cles, the electrons will suffer collisions and their average 
velocity will therefore decrease. This affects the space charge 
field and the current density. In the limiting case, where the 
anode-cathode separation is much larger than a mean free path, the 
velocity of the electrons is limited by their mobility in the gas. 
Cobine [55] derives an expression for the current density in a 
parallel-plate diode with a mobility-limited current: 

j = 9.92 X 10 (amp-cm (A.l) 

where d is the electrode separation in cm and is the mobility. 
The derivation assumed the electron velocity is proportional to 
the field, i.e. the mobility is field independent. This result 
means that the current is proportional to the square of the anode- 
cathode potential. A derivation for a cylindrical situation under 
the same assumptions yields the same dependence on plate voltage. 
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Although a general treatment has not been attempted, the similarity 

of the derivation to that of the Languir -Childs law suggests that 

the result is geometry independent. VVith the additional assumption 

that the mobility u is inversely proportional to pressure, the 

current in a mobility limited diode is given by 

2 

I a — (A. 2) 

P 



For weak fields, the mobility is not f ield- independent 
L 40 j and the drift velocity is given by 

v^a(S) . (A.3) 

In such a case, it is easy to show that the diode current for 
parallel plate and cylindrical diodes is given by: 






(A. 4) 



The field-dependent, the field- independent and any 
intermediate case as well as the transition from space charge to 
mobility limited operation are included in the general form 



I = Cp 



-''i'' “2 



(A. 5) 



where C, k^ and k^ are dependent on the situation and must be 
empir ica 1 ly determined . 

b. Sensitivity 

To determine the sensitivity of the diode as a micro- 
phone, consider the simple circuit shown in Figure 21. If the 
effect of the sound on the diode is observed by measuring the 
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Figure 21, Diode Circuit. A simple circuit which can be 
used to study sound with a mobility-limited 
thermionic diode. 




Figure 22, Diode Current in Dependence on Pres.su e ter a 
Typical Diode in Argon and Neon, 
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voltage across R, then the sensitivity, S, of the diode circuit 
can be defined as: 



6V 



S = - 



out 



6p 



(A. 6) 



Using equation (A. 5), it is easy to shov/ that this is given by : 



S 



V 



DC 



V 



1 + k 



DC 



2 V-V 



DC 



(A. 7) 



where V is the d.c. component of the output voltage. This 
expression shows that the sensitivity increases with the d.c. 
voltage across R, in an approximately linear relation, until the 
output voltage becomes comparable with the supply voltage V. The 
signal voltage V produced by an acoustic pressure wave p is there- 
fore given by 

V = ^. (A. 8) 

c. Frequency Response 

Equation (A. 7) is not correct for short wavelengths. 

Implicit in the discussion thus far has been the assumption that 

the finite size of the diode would not affect results. Actually 

the diode averages the pressure over its length. When fringing 

effects are ignored, the diode current is proportional to its 

length. Therefore a diode of length contributes a current 

6i given by . 

6i = I ^ (A. 9) 

where I is the current for a diode of length f , considered uniform. 
The total current consists of the d.c. component plus the a.c. 
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component caused by a pressure disturbance. The current fluct- 
uation is given by 

I = I p . (A. 10) 

^ ikx-iCOt 

For a pressure ivave given by p = > which is propagating 

in the axial direction of the diode, the total current I is given 
by integrating equation (A. 9) over the length of the diode, wherein 
the real part of the pressure must be used. 



- f 

• R Jl J 



COS kx dx. 



ther efor e : 



7' - S 

' - R Po 



o 

s in kf 
kl 



(A. 11) 



The frequency response of the di>de is therefore 



d. Tempei u t ur e Effect 



The equations developed thus far indicate that the 
diodes can be calibrated in terms of absolute pressures by extra- 
polation from the zero frequency case, i.e. the d.c. characteristic 
of the diode. This requires that the pressure be varied and the 
current measured as a function of plate voltage at the pressures 
used. The pressure measurement is usually made at some point 
removed from the diode. Since the filament is hot, the effects 
of the thermal gradient between the diode and the measuring point 
must be considered. 

Many properties of particles in gases which are fre- 
quently expressed as functions of pressure are in fact primarily 
dependent on the particle density. In such cases, temperature 
effects can be compensated for by correcting the pressure /or 
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temperature. For example, mobility is given as a function of 
pressure even though it is a function of density. 

n = |J.(p) 



If the data is for pressures at temperature and 

of the point of concern is > smd the pressure is 

effective pressure, p^, for entry into the data is 

T 



P 



e 




the 

P2’ 



temperature 
then the 



Now consider the effect of 
on the diode characteristics, starting 



the temperature difference 
from Gauss’s law: 



and the expression 



for 




o 

current density 



(A. 12) 



a weak field so that 



gas. Consider cylindrical 
symmetry and negligible end effects. Since the current I is a 
constant, and related to j by geometry, the density can be expressed 
in terms of the current and other factors then substituted into 
Eq. (A. 12). The result can be written in cylindrical coordinates 
as : 




j = pv. 

Assume electrons moving in a gas under 




where is a constant peculiar to the 
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This can be re-expressed as; 






^ A r — I = — r 

dr / dr dr / 2n£ 



(A. 13) 



The pressure p can be expressed in terms of the 
pressure, p^^, at the measuring point as 

p = p/(T(r)) = p^/(r) . 

With this substitution, Eq. (A.l6) can be integrated 



y di / 2 2n£ 3 ^ ^ 



Therefore 



^ = i 3 _i_ ^ F(r ) 

dr ] 2 2 nX c. 



2/3 



wher e 



F(r ) = ( / 2(r ' )r 'dr 

O 



from which the result 



' 3/2 



dr 



can easily be obtained. This can be written as 

tr3/2 



I = A p . 

^ "2 



A review of the steps taken to arrive at this result shows that 
even had some other dependence on pressure been used, the functional 
dependence does not change because of a temperature effect. 
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A theoretical treatment has given several predictions 



as to the operation of a mobility-limited thermionic diode as a 
microphone. The results may be summarized as follows. 

1. The diode d.c. current will vary with pressure and plate 

voltage as , 

“^1 2 
I = Cp V 

where % ^ ^ 1 and ~ ^ ^ 2 . 

2. For a simple circuit like Figure 21, the sensitivity, S, 
is given by 



V 



S = 



DC 



V 



1+ k 



DC 



2 V-V 



DC 



3 - 



sensit iv: 



m i c r op h on c w ill fa 11 off V7 i t h v/a v e — 



length as 



sin (2nf/X) 
(2Tif /A ) 



2. Experiment 

An experiment was performed using several diodes which 
investigated the first two predictions. The frequency response 
was not tested. 

a. Experimental Arrangement and Procedure 

The diodes used were of similar construction. They 
consisted of a tantalum cylinder, 1 cm long by 1 cm diameter, as 
plate and a tungsten or thoriated tungsten filament as cathode. 

The diodes were mounted several decimeters from the tube ends and 
the sound source, in a tube 4 cm in diameter. A current-modulated 



glow discharge acted as sound source. The a.c. signal was 
measured using a circuit like Figure 6 with a differential amp- 
lifier to amplify the signal which was then observed with an 
oscilloscope or with a boxcar integrator . 

After outgassing the system and activating the fila- 
ments, space charge limited operation in a high vacuum was verified. 
This was done by measuring the plate current as a function of plate 
voltage to a value of 3 ^ volts, which was twice the highest ex- 
pected operating voltage. If a log I, log V plot gave a straight 
line, the selected operating point was considered valid. If not 
and the current dropped below a straight line at 30 volts, then 
the diode was considered emission limited and appropriate steps 
taken. If the operating point was valid, then that filament 



voltage was maintained throughout the calibration. For the re- 
mainder of the calibration, the pressure was changed and the I 
versus V measurement made. For this portion of the calibration, 
if the log log plot showed that the current at higher voltages v;as 
above the linear fit to the lower values, than ionization was 
taking place and subsequent operation below these values of plate 
voltage was required. 

After calibration, the a.c. response to various pres- 
sure waveforms was noted. The sensitivity was measured with a 
boxcar integrator set to read the same point in the signal as the 
d.c. signal was varied by changing the resistance of the load 
resistor. The sound source amplitude v/as kept constant during the 
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measur ements . 



b. Noise 



The diode was found to be susceptible to electrical 
noise from external and internal sources. The external noise v/as 
reduced to an insignificant level by shielding with heavy copper 
foil. The foil was wrapped around the tube at its surface, ex- 
tending beyond the diode approximately 5 cm at both ends. The 
shield was grounded with a separate head. All leads to the diode 
were also shielded. 

In addition to the external noise, the diode v;as sub- 
ject to an internal noise which seemed due to spurious ionization. 
This noise v/as characterized by fast rising positive spikes, re- 
presenting a sudden increase in current, followed by a more gradual 
fall. These spikes always appeared when contamination was suspected, 
c. D.C. Characteristics 

Dayton et.al. L 53 l reported results for a diode oper- 
ating in whereas Argon and Neon were used in this experiment. 

The results obtained were similar to those of Dayton. A typical 
plot for Argon and Neon is shown in Figure 22 . The values of k^ 
and k^ obtained are displayed in Table IV along with information on 
N^ taken from Dayton^s paper. The values shown are the averages 
computed from data for several diodes for the pressure range 5 to 
30 tor r . 

TABLE IV 



Gas 






^2 


“3 


Ne 


.60 + .01 


1.67 


+ .04 


1.5 


Ar 


.63 + .02 


1.03 


+ .09 


1 


^2 


.62 


1.98 




1.7 
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d. Sensitivity 

A plot of signal voltage versus d.c. voltage for con- 
stant acoustic pressure constitutes a measurement of the sensitivity 
in arbitrary units. Figure 23 shows the results of such a plot. 

The solid curve was calculated using equation (A. 17) with and 
the amplitude determined by a least squares fit to the data points. 
The upper half of the points were neglected in the fit because of 
a definite shift in the slope of a linear plot of the data when the 
plate voltage became comparable to the filament voltage. 

3 . Discussion of Results 

The value of is very similar in all three gases but 
there is a large difference in the values for k^ • Since the 
theoretical treatment assumes that the pressure dependence arises 
from the similarity parameter E/p, if k is the same, then k 

JL ^ 

should be the same. As this is not the case, and the analysis 
has also shown that the temperature gradient should not affect 
this, then there must be a pressure dependence other than E/p. 

Von Engle points out that molecular gases , such as 

N^j have large in-elastic collision cross sections which causes 
their mobility to be constant and their drift velocity to vary 
linearly with E/p. This is related to the manner in which If, the 
average fraction of electron energy lost in a collision, increases 
with E/p. Until inelastic collisions start to take place, re- 
mains constant, afterwards it increases with E/p. The relation 
between the dependence of 5f on E/p and the dependence of the 



1^3 



9 



8 




Figure 23. . Diode Sensitivity in Dependence of d.c. 

Signal. Argon at pressure of 7*35 torr . 
Vertical scale is arbitrary. The ordinate 
is the d.c. component of the signal voltage. 
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mobility-limited diode current on plate potential is a complicated 

one and would be difficult to establish quantitatively. However 

there is a correlation between the experimentally determined values 

of and the slope of the log If - log E/p curves as figure 63 of 

Von Engle [4^1 • For simplicity, k is defined by 

3 




which is used merely for comparison. The values of k^ determined 
from Von Engle axe shown in Table iv^or comparison. 

4- Conclusion 

The mobility-limited diode nicrophone should serve to fill 
a need in plasma-acoustic studies. It is simple and inexpensive 
to construct and requires no complicated circuitry. By virture of 
changing the d.c. signal, sensitivities of different diodes can be 
matched for direct comparison of signals. It is sensitive. 
Sensitivities higher than .27 mV/microbar (3^0 mV/torr ) at 5 torr 
ambient pressure, can easily be attained. It is transparent to 
sound and offers little interference to light in the axial direction. 
Calibration at zero frequency is simple and easy to perform. 

It is susceptible to noise from contamination. This is a 
drawback which conventional microphones do not have. Its use 
would be limited to inert gases or those where chemical activity 
with the hot tungsten filament would not be expected. Even with 
these drawbacks, a diode microphone can be used for many experi- 
ments where other techniques can not. 
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APPENDIX B 



THE USE OF BUTT SEALS FOR FEEDTHROUGHS 

In order to seal electrode supports and feedthroughs into 
glass with a minimum disturbance of the glass surface, a special 
procedure was required. Kovar wire, 70 mil. diameter, was used 
for all supports and feedthroughs. Tantalum sheet 5 mil., was 
used for the electrodes. Diode filament leads were wrapped with 
fine rhenium wire for ease in spot welding. The electrode supports 
went along the electrode in the axial direction and were spot 
welded in several places. In addition, a small cut in the elec- 
trode was made at the edge and the support wire was bent back and 
crimped to hold the electrode in a mechanical grip. The small 
cut was made to prevent the wire from extending beyond the edge of 
the electrode. The materials w^ere cleaned prior to assembly and 
the components were cleaned again gust before sealing in accordance 
with procedures recommended by Rosebury C 56 I. 

For final assembly and sealing, the components of an element 
were assembled on a special -jig. The jig was designed to hold the 
support wires and feedthroughs in a plane which would be normal to 
the axis of the glass tubing. After assembly of the components the 
jig was attached to one piece of the glass tubing, earlier mounted 
in the lathe. The jig was placed so that the feedthroughs were 
flush against the edge as shown in Figure 24. The edge used for 
the seals worked best if it was diamond cut. Fire cut edges were 
susceptible to leaks. The second piece of tubing w^as matched with 



the edge of the first as shown in Figure 25 . A low pressure 
purge of nitrogen was used to prevent the flame from penetrating 
the tube and contaminating the surfaces. After the purge was on 
and the two pieces and jib were turning in the lathe, heat was 
applied and the two ends pushed toward each other. When suffi- 
ciently hot, the glass v;ould flow around the leads and the two 
ends would butt together, forming a butt seal. A small ring of 
7052 was then beaded around the base of the leads, for insurance. 
The complete tube was annealed in an oven following the recommended 
temperature cycle. For the tube used in the experiment, due to the 
nature of the main discharge filaments, annealing was done with 
the tube being evacuated by a diffusion pump. Enuring pre- 
liminary work, scvera .1 sea.ls were not annealed in the oven, but 
merely flame annealed on the lathe with no later cracking. 
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Figure 24. Diode in Process of Assembly, 
Showing Placement in Jig. 




Figure 25. Diode in Process of Assembly, 
Showing Assembly in Lathe* 
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APPENDIX C 



TAPE RECORDER SPECIFICATIONS 



FM Recor d/Repr oduce 
Input Sensitivity: (High Level) 



Input Impedance: 



(Low Level) 
(High Level) 
(Low Level) 



1 to 20 volts peak-to-peak 
for j+i 40 % deviation 

0.030 to 1 volt peak-to-peak 
for ^ 40 % deviation 

20k ohms or greater , unbalanced 
to ground 

50 k ohms or greater, unbalanced 
to ground 



Carrier Frequency, Response, Signal-to-Noise Ratio: 



High Level 


Center Frequency 


Response 


S/N 




.375 ips 


c 

0 


hz 


DC to 100 hz + 1 db 


35 


db 


3.75 ips 


5 


khz 


DC to 1 khz 1 db 


38 


db 


37 • 5 ips 


50 


khz 


DC to 10 khz + 1 db 


42 


db 


Low Level 












.375 ips 


0 

0 


hz 


5 to 100 hz + 1 db 


35 


db 


3.75 ips 


5 


khz 


5 hz to 1 khz + 1 db 


38 


db 


37 . 5 ips 


50 


khz 


5 hz to 10 khz + 1 db 42 


db 


Gain Drift: 






2 % or less of fall 
scale 


( t emp . 

35° to 



Center Frequency Drift 

Linear ity : 

Output Level: 

Output Impedance: 
Harmonic . Distortion ; 



2% or less of full 
sea le 



line voltage 
range 110 
vac to 125 
vac ) 



1 % of best straight line through 
zer o 

2.8 V p-p for ^ 40 % deviation 
into 1,000 ohms 

50 ohms or less 

1 . 5 % total (signal 1 v p-p, 

3,000 hz, at 37.5 ips) 



149 



Output Filter 



Cutoff Frequency; 

Fr equency Response : 



Overshoot : 



Transport 



Tape Speeds: 



100 hz , .1 khz, 10 khz, 100 khz 

Filter set for flat frequency 
response, ^ 1 db -- Filter set for 
flat phase response -4 to -8 db 
at cutoff. 

Filter set for flat frequency 
response, 15% raax. -- Filter set 
for flat phase response, 2% max. 



37.5, 3.75, 0.375 ips 



Speed Accuracy: 



+ 0.2%, long-term 
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COMPUTER PROGRAMS 



I. DATA REDUCTION 

AFTER THE DATA WERE STORED ON TAPE IN 32-BIT 
BINARY FORM, FIVE PROGRAMS WERE USED SEQUENTIALLY TO 
PROCESS THE DATA. DISK STORAGE WAS USED FOR INPUT AND 
OUTPUT EXCEPT THE FIRST PROGRAM USED TAPE FOR INPUT DATA 
EACH PROGRAM MAY BE THOUGHT OF AS A STEP IN THE PROCESSING 

IcSTEP ONE — AVERAGE — 

C THIS PROGRAM TAKES THE RAW DATA AND BEGINS THE 

C PROCESSING. the DATA HAS BEEN PRESCREENED AND INFORMATION 
C STORED IN AN ARRAY INDICATING WHICH SETS OF WAVEFORMS ARE 
C ACCEPABLE FOR AVERAGING. THIS ARRAY IS CALLED *KF». ALL 
C ACCEPTABLE WAVEFORMS ARE AVERAGED AND THE STANDARD DEVIA- 
C TION OF EACH DETERMINED FOR EACH DIFFERENT EXPERIMENTAL 
C RUN. THIS INFORMATION IS STORED FOR USE BY THE NEXT 
C PROGRAM. 

C 

c 

DIMENSION DATA(3600) ,AVG(3600) ,KSTAR(12) ,KF(100,12) 
DIMENSION STDEV( 360Q),DEV(3) 

REWIND 4 
REWIND 2 

DATA KF,NF,KSTAR/1200«0, 1,3=!‘2,3, 1,1, 3, 1, 1, 2/ 

C 

READ( 3)KF 

5 NR=0 
C 

KG=0 

2 DO 7 I=1,360C 
STDFvm=0. 

7 AVG(I)=0.0 
1 NR=NR+1 

READ(2,ERR=1,END=50) DATA 
IF(KF(MR,Nr).EQ.O) GO TO 1 
C 

DO 4 1=1,3600 

STDEVI I ) = STDEV( I ) + OATA( I >*>1=2 
4 AVG( I) = AVG( I )+DATA( I) 

KG = KG-*-l 
GO TO 1 
C 

50 DIV=KG 

DO 6 1=1,3600 

STDEV( I )=SQRT( ( STDEV( I )-AVG( I ) ^«2/ D I V ) / ( D I V- 1 . ) ) 

6 AVG(I )=AVG( I )/DIV 
C 

DO 8 J=l,3 
DEV(J)=0.0 
DO 9 1=1,1200 
I I=I+( J-l)-1200 
9 DEV( J)=DEV( J )+STDEV( II )v*2 

8 DEV( J)=SQRT(DEV( J) /1200o ) 

C 

WRITE(4) AVG,DEV 

WRITE (6, 200) NF 

WRITE (6, 203) DEV 

WRITE! 6, 201) AVG 

WRITE(6,202) ( K F ( 1 , NF ) , I =1 , NR ) 

NF=NF+1 

IF(NF.LE.IO) GO TO 5 
C 

WRITE(13) KF 
C 

2C0 FORMAT! '1 FILE NO. *,I3) 

201 F0RMAT(60! 1 X , 1 0 F 12 « 4 , / ) ,/) 

202 FORMAT! ' 0* , 10! lOX, 101 5 , / ) ) 
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no 



203 FORMAT {' 
STOP 
END 



AVERAGE DEVIATIONS 



• ,3F10. 5) 



2. STEP TWO — SCREEN AND RE-AVERAGE — 

C THIS PROGRAM USES THE AVERAGE WAVEFORMS AND THE 

C STANDARD DEVIATION CALCULATED IN THE PREVIOUS PROGRAM TO 
C SCREEN OUT DATA WHICH LIES OUTSIDE OF TWO STANDARD DEVIA- 
C TIONS OF THE AVERAGE. IT RECOMPUTES THE AVERAGE AND 
C STANDARD DEVIATION WITH THE SCREENED DATA. 

C 

C 

DIMENSION DATA (3 600) ,AVG(3600) ,KSTAR(12) ,KF(10C,12) , 
>!'KG(36C0) ,DEV(3) , STDEV(3600) 

REWIND 2 
REWIND 4 
NF=1 

READ(3) KF 
REWIND 3 
STD=2. 

C 

5 NR^O 

READ! 4) AVGfDEV 
DO 2 1=1,3600 
STDFV( I )=AVG( I )^f^2 

2 KG( I )=1 
1 NR=NR+1 

READ(2,FRR=1 ,EM0=50) DATA 
IF(KF( NR,NF ).EQ.OO) GO TO 1 
C 

KN=0 

DO 3 I=1,36C0 
j=T 200 + 1 
01 V=KG ( I J 
COM=AVG( I )/DIV 
TEST=ABS( COM-DATA( I ) ) 

TEL=STD^DEV( J) 

IFITESTcGT.TEL ) GO TO 3 
STOEV( I ) =STDEV ( I )+DATA ( I )*^2 
AVG( I )=AVG( I )+OATA( I > 

KN=KN+1 

KG( I)=KG( I )+l 

3 CONTINUE 
K>=(NR,NF) = KN/36 
GO TO 1 



50 DO 6 I=1,36CC 
OIV=KG( I ) 

STDEV( 1 )=( STOEV( I)-AVG( I ) 4=^2/ D I V ) / { D I V- 1 . ) 
IF( STDEV( I ) .LT.0.0) WRITE(6,203) I,STDEV(I) 
STDEVI I )=SQRT( ABS( STDEV( I ) ) ) 

6 AVG( I ) =AVG( I )/Dl V 

DO 8 J=l,3 
DEV( J) =0.0 
SUM=0« 

DO 9 1=1,1200 
I I = I + ( J-1)*120C 
SUM = SUM+KG( 1 I ) 

9 DEV{ J)=OEV( J)+STDEV( II )*«2*KG( I I ) 

8 DEV( J)=SQRT(DEV( J)/SUM) 

WRITE(9) AVG.DEV 
WRITE(6,2C0) DEV 
WRITE (6, 201) AVG 
V/RITE ( 6, 202 ) ( K F ( I , NF ) , I = 1 , NR ) 

NF=MF+1 

IF(NF.LE.IO) GO TO 5 
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WRITE(13) KF 

200 FORMAK’l DEVIATIONS ARE S3E12.3) 

201 F0RMAT(60{1X,10F12.4,/),/) 

202 FORMAT! 'O' , 10( ICX, 1015,/) ) 

203 FORMAT!' STOEV SQUARED IS NEGAT I VE ' , 1 5 , E 1 2. 2 ) 
STOP 

END 



3. STEP THREE — DELAY , OVERLAY AND ZERO — 

THIS PROGRAM TAKES THE AVERAGED WAVEFORMS AND 
DELAYS THEM, REDUCING THE NUMBER OF POINTS FROM 1200 TO 
600 IN EACHo IT ALSO ZEROES THE BASE AND OVERLAYS A 
FALSE SIGNAL DUE TO AN ELECTROSTATIC INTERACTION. THE 
RESULTS ARE STORED FOR THE NEXT STEP, WHICH IS THE FOURIER 
TRANSFORMATION. 



DIMENSION DA! 1200, 3) , I ND ! 3 , 10 ) , S I G ! 600 , 3 ) ,DEV!3) , AVG!3 
1 IPOINT! 2 ,3, 1C) 

EQUIVALENCE ! D A ! 1 , 1 ) , S IG ! 1 , 1 ) ) 

DATA I PCI NT/0, 600,0,600,0,600, 2AC , 60C , 1 2 0 , 2 70 , 0 , 288 , 24 
*96,288,0,288,240,600, 121,238,0,288,240,600,204,276,0,2 
*240,600 ,216,288, 0,264,240,600, 108,275,0,250,240,600,10 
*0 , 24C , 2 40 , 6 00 , 1 20 , 2 20 , 0 , 240 , 240 , 6C 0 , 1 20 , 2 40 , 0 , 2 40 / 

KF=0 
MF = 0 

REWIND 4 

READ! 5, 100) ! ! I ND ! J , I ) , J= 1 , 3 ) , I =1 , 10 ) 

301 NF=NF+1 

READ!4,END=305) DA, DEV 
KF=KF+1 

IF!NFoEQoD GO TO 312 
A=URN!0 ) 

DO 311 1=100,240 
JI=130.*URM!l)+240 

311 DA! I, 2) =DA! JI, 2) 

312 CONTINUE 

DO 304 J=l,3 
KC=IND! J,KF ) 

DO 302 1=1,600 

302 SIG! I , J )=OA!KC+I , J ) 

304 CONTINUE 

DO 321 J=1 , 3 
AV6! J)=0.0 
IP1 = IP0INT! 1 ,J,NF)- 
IP2=IP0INT! 2, J , NF) 

DO 322 I=IP1,IP2 

322 AVG!J)=AVG!J)+SIG!I,J) 

DIV=IP2-IP1 

321 AVG! J)=AVG! J)/DIV 

no 323 J=l,3 
DO 324 1=1,600 

324 SIG!I,J)=SIG! I,J)-AVG!J) 

WRITE 16,200) KF,J 
00 303 1=1,589,12 

303 WRITE!6,202) ! SI G ! I +KM- 1 , J ) , KM=1 , 1 2 ) 

323 CONTINUE 

WRITE!8) SIG ,DEV 
GO TO 301 



305 STOP 

100 FORMAT! 1215) 
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200 FORMAK'l R.TRIPP FILE ',13,’ CHANNEL', lA) 

201 FORMAT! 20A4) 

202 FORMAT! 1X,12F10.4) 

END 



4. STEP FOUR — TRANSFORM — 

THIS PRCGRAM DETERMINES THE SPECTRAL DECOMPOSITION 
OF THE WAVEFORMS AND THEIR STANDARD DEVIATIONS. GAIN AND 
DIODE SENSITIVITY ARE CORRECTED FOR. 



INITIALIZE VALUES 

DIMENSION SIG! 400,3 ) , SPECT! 600,3) ,T!600 ) , FREQ(60 ) , 

1 SHO! 200) , BI AS!3) ,GAIN! 3) , SENS! 3) , IP!3) , 

2 SDEV! 40, 3 ) ,SREQ!60) , DEV!3) , INV! 2CC) ,S!20G) 

• EQUIVALEfJCE ! S I G , SPECT ),! SHO , SDEV ) 

COMMON SIG, SDEV, S, INV 

DATA PI , GA I N, N PNTS, NF , PERIOD/3 o 1415 92 65 , 3=«= 50. ,600,0, 

1 0.00512/ 

DO 300 I=1,NPNTS 

300 T! I )=I 

DO 301 1=1, 40 
FREQ! I) =!I-i )/PERIOD 

301 SREQ! I )=SQRT(FREQ! I ) ) 

READ DATA FOR NEW DATA SET 

17 CONTINUE 
NF=NF+1 

READ!8, END=9o) SIG, DEV 
RFAD!5,101) SENS 
WRITE!6,101) SENS 

CORRECT FOR AMPLIFIER GAIN 
AND DIODE SENSITIVITY 

DO 705 J=l,3 
FACT=GA IN! J )*SENS! J ) 

DEV! J)=DEV! J)/FACT 
DO 303 I=1,NPNTS 

303 SIG! I,J)=SIG! I ,J)/FACT 

DISPLAY WAVEFORM ON PRINTER 

DO 305 1=1,199 

305 SHO! I ) = !SIG!3>!^I-2, J)+SIG!3=!‘I,J)+SIG!3’i=I + 2,J) )/3. 
SH0{2f)C) = SlG!6C0, J) 

WRITE!6,200) 

705 CALL PLOTP! T, SHO, 200,0) 

IF THIS IS THE FIRST DATA SET, DETERMINE CORRECT 
GAIN FOR REMAINDER OF DATA SETS. 

IF!NO,EQ.l) GO TO 313 

DETERMINE SPECTRUM AND STD. DEVIATION 

CALL SPCTRA 
DO 304 J=l,3 

304 SOEVII, J) = SQRT!SDEV! 1, J)=<^«2 + 0EV!J)^''’^2) 

18 WRITCIlO)! !SPECT!I ,J),SDEV!I,J) ,J=1,3 ) , FREQ ! I ) , I =1 , 25 

DISPLAY RESULTS AND PASS THEM TO NEXT 
PROGRAM 

DO 307 J=l,3 
WRITE! 6,200) 
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DO 309 1=1,40 

SIG( 1 + 40, J)=SIG( I, J)-S0EV( t , J) 
IF(SPECT(I+4G, J ).LT.0.0) SPECK I +40, J) =0,0 
SIG(I+B0, J)=SIG( I, J)+SDEV( I, J) 

309 CONTINUE 

DO 308 KJ=1,3 
KM=121 -K 

CALL F'L0TP(fREQ,SIG(KM,J) ,40,KJ) 

308 CONTINUE 
307 CONTINUE 



1=0 

DO 302 K=l,27 
1 = 1+1 

WRITE! 6,210) K, (SPECT(K,J) ,S0EV(K,J) , J=l, 3),FREQ(K) ,SR 
RATI0=SPECT (K, 3 ) /SPECT ( K,2) 

SPECT( I , 2) = AL0G( SPECT(K, 2) ) 

SPECT n ,3)=RATI 0 
■ SPECT ( I , 1 )=AL0G( RATIO) 

302 CONTINUE 



WRITE (6, 200) 

CALL PL0TP(FREQ,S'^F.CT(1,3) ,27,0) 

CALL PL0TP( SREQ, SPECT (1,1) ,27,0) 

GO TO 17 

CALIBRATE GAIN 

313 CALL RHARM(SPECT(1,2),8,INV,S,IER) 

SPECT( 1 ,?)= SQRT( SPECT( 53,2) 2 + SPECT (54,2) =f'>;-2 ) 

GAIN(2 ) = GAIN( 2 ) ^J'SPECT ( 1, 2) 

CALL RHAFM( SPECT (1 ,3) , 8 , 1N\/ , S , I ER ) 

SPECT (1, 3) = SQRT( SP ECT ( 53 , 3 ) =^‘^2 + SPEC T ( 54 , 3 ) *« 2 ) 

GA!N(3 ) =OA T N(3 ) ^SPFCT ( 1 , 3 ) 

WR I TE ( 6 ♦ 10 1 ) G A I N , S DE C T ( 1 , 2 ) r SP ECT ( 1 , 3 ) 

GO TO 17 



FORMAT STATEMENTS 

100 FORMAT! 20A4) 

101 FORMAT (3F1C. 5) 

200 FORMAT (IHl) 

202 FORMAT(lHl,/,‘ TIME ANO SIGNAL • , / , 50 ( 1 2F 10, 4 , / ) ) 

203 FORMAT! IHl,/, TIC, • SPECTRAL ANALYSIS OF SIGNAL FOR ', 

I'OIOOE • ,12,/, T6, ' INDEX FREQUENCY COSINE •, 

2'SINE AMPLITUDE PHASE') 

204 FORMAT! 1 10 , 1P3E 1 1. 3 , 0P2F 1 1 . 4) 

206 FORMAT! • UVL=',E10.2) 

207 FORMAT!’ FILE KI4) 

210 FORMAT! I 5 , 8E 1 1 . 3 , I 5 ) 

211 FOR'^AT! 20A4) 

500 FORMAT !1H ,T10,'LESS THAN TEST') 

99 STOP 
END 



SUBROUTINE SPCTRA 

THIS SUBROUTINE DETERMINES THE FOURIER TRANSFORMS 
OP THE WAVEFORMS FOR FIVE START POINTS. K THEN AVERAGES 
THEM AND CALCULATES THE STANDARD DEVIATION OF EACH ONE. 



SUBROUTINE SPCTRA 
IMPLICIT C0MPLEX*8!0) 

DIMENSION SPEC ! 51 6 , 3 ) , S ( 200 ) , I NV( 200 ) , M ! 3 ) , A ! 2 , 2 58 , 3 ) 
1 SIG!600,3),S0EV!4D,3) ,AVG!40,3) 

COMMON S I G , SDE V , S , I NV 
EQUIVALENCE (SPEC, A) 
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DO 201 1=1, ^0 
DO 202 J=l,3 
AVG( I , J )=0.0 
202 SDEV( I , J)=0.0 
201 CONTINUE 
C 

c 

DO 706 KT=1,5 
NL=10«( KT-l ) 

DO 705 J=l,3 
DO 7C7 1=1,512 
SPEC( I , J ) = SIG( NL + I , J) 

7C7 CONTINUE 

CALL RHAPfM SPEC( 1,J),8,INV,S,IER) 
IF(IER.NE.O) WRITE(6,10C) lER 

705 CONTINUE 
C 

C 

DO 302 J=l,3 
DO 301 1=1,40 

RSPEC=A( 1, I , J) *=f'-2+A(2, I , J ) «*2 
SDEV( I, J)=SP5V( I,J)+RSPEC 

301 AVG(I ,J)=AVG(1,J)+S0RT(RSPEC) 

302 CONTINUE 

706 CONTINUE 

DO 304 J=l,3 
DO 303 1=1,40 

SDEV( I, J)= ABS { { SDEV( I , J )~AVG( I , J ) **2=!'% 2 ) ) 
SDEV( I , J) = SOPT{ SDEV( I , J) )>■''. 2 5 

303 SIG( I , J )=AVG( I , J ) + .2 

304 CONTINUE 
RETURN 

100 FORMAT (• HARM ERROR IERRT=',I3) 

FMO 



5. STEP FIVE — TRANSFER FUNCTIONS— 

C THIS PROGRAM TAKES THE FIRST 25 COMPONENTS OF THE 

C SPECTRAL DECOMPOSITION OF THE WAVEFORMS AND USES THEM TO 
C DETAIN THE ACOUSTIC TRANSFER FUNCTION. IT ALSO DETERMINES 
C THE SOURCE TRANSFER FUNCTION AND DISPLAYS THE DATA ALONG 
C WITH A CURVE REPRESENTING TH'^ PREDICTED FORM , V.'HICH WAS 
C FITTED TO PASS THROUGH ONE POINT. THE STANDARD DEVIATIONS 
C OF THE RESULTS ARE ALSO CALCULATED. 

C 

IMPLICIT REAL^i'B ( D ) , COMPLEX’i'8 (Q) 

REALMS TITLE(10)/7*’ ', ’TRIPP ’,2^' •/ 

DIMENSION DX{ 25) ,DY{ 25) ,DF{ 25 J , DW{25) ,DELY{25) , 

1 DB{20) ,DSB(20) ,X(3) , SPECT(25,3) ,FREQ(25), 

2 XX (2) , YY (2),SDEV{ 25,3 ),Z( 25,3 ),UL(2), 

3 SREQ{ 80) , YZ ( 30) ,USRCE( 25) 

C 

DATA DW,X/ 254=1. DO, 0.0, 98. 7, 251. 3/, 

1 XX, YY/C. 0,4500. ,10. ,0.0/ 

K=25 

NF=0 

10 CONTINUE 
NF=NF+1 

READ! END = 99) { ( S PECT ( I , J ) , SOE V { I , J ) , J = 1 , 3 ) , 

1 FREQI I ), 1=1,25) 

DO 300 1=1,25 

300 USRCEI n=SPECT( I,3)/SPECT(I ,2) 

WRITE(6,2C3) (USRCEI I ) ,FREQ( I ), 1=1,25) 

C 

C CONVERT TO DOUBLE PRECISION 

C 

L=0 

DSUM=0.0 
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DO 301 I=2tK 

FSPQ=(SDEV(I,2)/SPECT(I,2) )«*2 
1 +(SDEV( I,3)/SPECT(I,3) )«*2 
L = L + 1 

DY(L)=ALCG(USRCE( I) ) 

DW( L)=1.D0/FSP0 
OSUM=DSUM+DW(L ) 
nX(L)=SQRT(FREQ( I ) ) 

301 CONTINUE 

IF(L.LT.4) GO TO 10 
DL=L»^50 
DO 311 I=ljL 
311 DW( I )=DW( I i^DL/OSUM 
C 

C DETERMINE ALPHA 

C 

CALL LSQPL2( L, 1 » DX , DY , OW , DF » DEL Y , DB t DS B , T ITLE ) 
ALPHA=-SNGL (DB(2))/(X(3)-X(2)) 

UALPHA=SMGL(DSB(2) ) / ( X ( 3 ) -X ( 2 ) ) 

WRITF(6,205) ALPHA, UALPHA 
C 

C DISPLAY LEAST SQUARES FIT 

C DISPLAY TUBE TRANSFER FUNCTION 

C 

ZINC=FREQ( K)/79. 

FA=DB( 2 ) 

YZn )=DB( 1 ) 

YZm = FXP(YZ(l ) ) 

. SRFQ( 1)=0. 

DO 304 1=2,80 
SREQI I ) = ( I-l )=f^ZINC 

304 YZ( I )=YZ ( 1 )-"EXP ( F ASSORT! SR EQ( I ) ) ) 

C 

CALL PLOTP(SREO,YZ, 80,1 ) 

CALL PLOTP(FREQ,USRCE,25,3) 

GOTO (11,10,10,10,10,11,10,11,10) ,NF 
1 1 rnwT] Ml 

C determine SOURCE AND SOURCE UNCERTAINTY 

C 

UL( 1)=0. 

UL(2)=G. 

WRITE (6, 204) 

Al=251.3/152.6 
A2=-98.7/1F2.6 
DO 303 I=1,K 

Z( I,1)=SPECT(I ,2)«^A1*SPECT(I ,3)*^ A 2 
USRCE ( I ) = ( A1 ’fSDtVI I ,2 ) / SPECT (1,2) )-'^*2 
1 +(A2*SDFV( I,3)/SPECT(I,3) )*’f^2 

USRCF( I)=2( I , 1 )^SQRT(USRCE( I) ) 

Z(I,2) = Z(I,n+USRCF(I) 

Z(I,3) = Z(I,1)-USRCE(I ) 

UL(2)=AMAX1(Z( I ,2) ,UL( 2) ) 

V1RITE( 6,201 ) I,Z(I,1),USRCE(I),FREQ( I) 

303 CONTINUE 
C 

C DISPLAY SOURCE SPECTRA 

C 

WRITE(6,202) 

CALL PL0TP(XX,UL,2, 1) 

CALL PLOTP( FPEO,Z( 1 ,1 ) ,K,2) 

CALL PLOTP (FREQ, Z( 1 ,2 ) ,K,2) 

CALL PLOTP( FREQ,Z( 1 ,3) ,K,3) 

C 

C DETERMINE TRANSFER FUNCTION AND 

C UNCERTAINTY 

C 

YY( 1 )=0.C 
DO 302 I =1 ,K 
TZ = Z( 1,1) 

Z( 1 , 1 ) = :/( I , 1 )/SPECT(I , 1 ) 

SDEV( I ,1 )=Z ( 1 , 1 ) -i^SQRT ( (USPCF( I )/TZ)=<'-2 
1 +(SDEV( 1, 1) /SPECT ( I ,1) )^-^0?) 
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Z( 1 ,2 )=Z( I , 1 )+SDEV( 1,1) 

Z(I ,3)=Z(I ,1)-S0EV( 1,1) 

WRITE ( 6, 201 ) I , Z(I , 1 ) , SDEV( I ,1 ) ,FREQ( I ) 
YY(l)r:AMAXl(YY(i),Z(I,l) ) 

302 CONTINUE 

YY(1)=YY(1)*1.2 

C 

C DISPLAY TRANSFER FUNCTION 

C 

WRITE( 6,200 ) 

CALL PLOTP( XX, YY,2 ,1 ) 

CALL PLOTPIFREQ, Z( 1 ,1 ) ,K,2) 

CALL PL0TP(FRE0,Z(1,2),K,2) 

CALL PLCTP( FREO, Z( 1 ,3) ,K,3) 

YS=0.75E-0A 

CALL DRAVv(K,FP':Q,Z(1,1) , 1 , 1 , T I TL E , T I TLE , 10 CO . , YS , 0 , 0 , 
1 2, 2, S, 6,0, LAST) 

CALL DRAW(K,FREQ,Z( 1,2) , 2 , 2 , T I TLE , T I TL E , 1 000 . , Y S , 0 , 0 , 
1 2, 2, 5^4,0, LAST) 

• CALL DRAW(K,FRcQ,Z( 1,3),2,2,TITLE,TITLE,1C00. ,YS,0,0, 
1 2, 2, 5, 4,0, LAST) 

C 

C DISPLAY PREDICTED CURVE FITTED TO DATA 

C 

READ(5,102) AA 
102 FORMAT! F 10. 5) 

FACT = STN( AA=}^1.0134E-C4^^FREQ( 14) )/FREQ( 14) 

FACT2 = Z( 14, D/FACT 
SREQ(l)=lo 
DO 306 T=l,80 

FACT=SIN(AA^1„C134E-04*SREQ( I))/SREQ( I) 

306 YZ( I )=FACT-FACT2 

CALL DRAW(80 ,SREQ, YZ, 2 ,G, TITLE, TITLE, 1000. ,YS, 0,0, 

1 2, 2, 5, 4,0, LAST) 

DO 305 I=1,K 
SREO( 1 )=FR?^0( I ) 

SRFO( 2 ) =FRFO( I ) 

SDEV(2,1 )=Z( 1,3) 

SDEV(1,1)=7(I,2) 

KK=I/l< + 2 

3C5 CALL 0RAW(2,SREQ(1 ) ,SDEV,KK,0, TITLE, TITLE, 100-). ,YS,0, 
1 0,2,2 ,5, 4,0, LAST) 

GO TO 10 
P9 STOP 

100 FORMAT! 20A4) 

101 FORMAT! 13) 

200 FORMAT! •! TRANSFER FUNCTION OF SOURCE') 

201 FORMAT! lOX, 13, 3E12.3) 

202 FORMAT! »1 SOURCE SPECTRA') 

203 FORMAT!'!^ RATIO FREQ', 

1 /,25! lOX, 2F1-2.5,/) ) 

204 FORMAT! 'O' ,T15, • SOURCE FREQ') 

205 F0RMAT!'0',' ALPHA I S ' , 1 Pi E 1C . 2 , ' +/ -' , ElO . 2 , / ) 

END 



/ ‘ 



158 



OOO OOO oooooooooooooo 



II. INGARn-SCHULZ ANALYSIS 



THE FIRST PART OF THE PROGRAM INITIALIZES VALUES. 
SYMBOLS FOR THE CONSTANTS OF THE GAS ARE: 

CD MOLECULAR DIAMETER 
CMN MASS OF ONE MOLECULE 
CG RATIO OF SPECIFIC HEATS 
VALUES OF THE INTEGRALS USED IN EVALUATING THE 
ELECTRON-NEUTRAL COLLISION FREQUENCY AND THE AVERAGE CROSS 
SECTION ARE CONTAINED AS A TABLE OF ENTRIES AS Y1 AND Y2 
RESPECTIVELY. EN IS THE VALUE OF ELECTRON TEMPERATURE 
USED AS ENTRY INTO THE TABLE. FOR THE COMPUTATIONS, THE 
HKSA SYSTEM OF UNITS ARE USED EXCEPT MASSES ARE IN GRAMS 
AND LENGTHS IN CENTIMETERS. 

REAL>t^4 X(5) ,XI (2) , STARPT(5) , L OWL I M ( 5 ) ,U PL IM ( 5 ) , XMN ( 5 ) / 
>J=NAME{ 6) 

DIMENSION XX(IOO) ,YY( 100) ,EE(40 ) , SIG(40) 

DIMENSION PQ{5) ,DS(5) 

CCMMON /1NNER/C(2C ) ,EN( ICC) , Yl{ 100) ,Y2( ICO) ,ARG( 10) , VA 

COMMON /RLOC1/GUT(2C)/BLOC2/MO 

PI=3. 1^159265 

CK=1.380622E-16 

CME=9.10955&E-28 

CA=.E2c>1771SE-8 

CQE=1 .6C21917E-19’ 

CF=8, 6542E-21 

EVALUATION OF CONSTANTS 

READ! 5, 101 ) CD,CMN,CG,CP ,NE,NAME 
DO 301 I =1,20 
OUT{ I )^0.0 
301 cm =0.0 

C( ] )=SORT(CG’:^CK/CMN) 

r ^ 21 = 1' /< 1 . Ai A2=!<p T«r.n«*2 ) 

C(3) = pi=i^o5*CA=t'=’)'2 

C(7)=CME/CMN 

C(9)=2.=«''CK/CME/PI 

C( A )=(CG-1 )>:^C( 7 )^*2*C(9) 5 

C( 5)=S0RT(8./P1/CG) 

C(6)=2.*C(3)=^^S0RT(C(9) ) 

C(8)=CG 
C( 10)=CR 

C( 11)=CQE/SQRT(CME-'^CE) 

C( 12)=f QE=:=A2/4./?I/CE/CK 
C( 13)=SQRT(CK/CME) 

C( 14) = SQRT( 1.667=<=CK/CME) 

C( 15) = 12.':^PI 
WRITE (6,201 ) C 

. ESTABLISH INTEGRAL TABLE 

READ! 5, 104) ( EN ( I ) , Y1 { I ) , Y2 ( I ) , 1= 1 , 100 ) 

C 

C AFTER INITIALIZATION, THE- MAIN PROGRAM CONTROLS 

C ACCESS TO THF SUBROUTINE 'IMGAPD'. THE SUBROUTINE 
C EVALUATES THE DISPERSION RELATION OF INGARD AND SCHULZ. 

C THE ARRAY 'GUT' CONTAINS I 

C THE ARRAY 'OUT' CONTAINS VALUES OF INTEREST CALCULATED IN 
C THE EVALUATION. THE CONTROL VARIABLE ‘KCC’ IS USED TO GIVE 
C A LARGE QUANTITY OF INFORMATION IN A SINGLE CALL TO INGARD 
C A VALUE OF KCC=1 CAUSES SEVERAL NAMELISTS TO RE PRINTED, 

C GIVING THE VALUES OF ALL INTERMEDIATE QUANTITIES. FOR ANY 
C OTHER VALUE, THE SECOND PART IS USED TO CALCULATE THE EXPE 
C ATTENUATION FACTOR FOR A SET OF PARAMETERS AND DETERMINE 
C THE UNCERTAINTY OF THE RESULT. 

C 

C 

KCC = 2 
J = 3 
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READ(5,103) STARPT 
CALL INGARD (STARPT) 

WRITE(6,209) 0UT(19) 

DETERMINE ATTENUATION CONSTANT FOR 
COMPARISON WITH EXPERIMENT 

C( 10=2. 5 
LR=1 

800 READ( 5, 103, END = 999 ) STARPT 
REA0(5,103) DS 
KCC=1 

DO 310 LQ=1,2 

IF(LO.GT.l) STARPT( 1)=6C0. 

WPITE(6,208) STARPT ,DS 
DO 304 MO=l‘3,17 
CALL INGARD(STARPT) 

KCC=2 

T0=0’JT(M0) 

DO 302 L=l,5 
TX = STARPT{ L ) 

STARPT (L)=TX*1,1 
CALL INGARD(STARPT) 

P0( L)=(OUT(MO)-TO) / (.1«TX) 

302 STARPT(L)=TX 
VAR=O.C 

DO 303 L=l,5 

303 VAR=VAR+(P0(L)*DS(L))**2 
VAR=SORT(VAR) 

WRITE(6,206) MO, TO, VAR, PO 

304 CONTINUE 
31C CONTINUE 

IF(LR.GT.l) C(10)=2. 

LR=LR+1 
GO TO 8 CO 



101 FORMAT(4E10.0, I2,6A1) 

102 FORMAT ( 2F10.0i 

103 FORMAT ( SElOoO) 

104 FORMAT (3E16.7) 

1C5 FORMAT! 31 5) 

201 FORMATdH , / , 4 ( IP 5E 12 . 3 , / ) ) 

202 FORMAT (IHl) 

203 FORMAT! IHl, /TIO, 'MINIMUM FOR OUT!', 12,') IS *,1PE9.1, 
=!=3X,' AT ',5E10.1) 

204 P0RMAT!1H0,T10, 'START POINT WAS ',IP5EIC.1) 

205 FORMAT! IHO, TIC, 'OUTPUT VALUES ARE : ' , / , 4 ! 1 P 5F1 2 . 3 , / ) ) 

206 FORMAT!' OUT!', 12,') IS • , E 12. 3, ' WITH STD DEV',F12.3, 

1 /,' PARTIALS ARE *,5E12.3) ' 

208 FORMATdH ,1P5E15.3,/) 

209 FORMAT!' SOURCE TERM ',1P1E12.3) 

214 FORMAT! 1M1,T3D, ' PLOT FOR ',6A1) 

<^99 STOP 
END 



SUBROUTINE INGARD!X) 



THIS SUBROUTINE SOLVES FOR THE REAL AND IMAGINARY 
COMPONENTS OF THE PROPAGATION CONSTANT *K' OF AN 
ACOUSTIC WAVE TRAVELING IN A CIRCULAR TUBE THROUGH AN 
IONIZED GAS, ACCORDING TO THE EQUATIONS DEVELOPED BY 
INGARD AND SCHULZ. 

VARIABLES AND THEIR SYMBOLS ARE*. 

ANGULAR FREQUENCY X!l) 

NEUTRAL DENSITY X!2) 



l60 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



c 



c 

c 



c 

c 

c 

c 

c 

c 



c 

c 

c 



ELECTRON DENSITY X(3) 

NEUTRAL TEMPERATURE X(4) 

ELECTRON TEMPERATURE X(5) 

VALUES OF CONSTANTS ARE COMPUTED IN THE MAIN 
PROGRAM AND PASSED IN COMMON, LIKEWISE THE PRODUCTS OF 
THE INTEGRATION SUBROUTINE ARE PASSED FROM THE MAIN 
PROGRAM IN COMMON, ALSO THE RESULTS OF THIS SUBROUTINE 
ARE PASSED BACK TO THE CALLING PROGRAM IN COMMON, THEY 
ARE STORED IN THE ARRAY 'OUT*. THE INDEX NUMBERS OF 
THIS ARRAY AND THE QUANTITIES STORED THERE ARE*. 

1 REAL PART OF N-SQUAREO 

2 IMAGINARY PART OF N-SQUARED 

3 REAL PART OF PROPAGATION CONSTANT 

4 IMAGINARY PART OF PROPAGATION CONSTANT 

5 MAGNITUDE OF ION VELOCITY 

6 PHASE OF ION VELOCITY 

7 MAGNITUDE OF ELECTRON VELOCITY 

8 PHASE OF ELECTRON VELOCITY 

9 RATIO OF OUT (4) TO THAT WITHOUT IONIZATION 

10 IMAGINARY COMPONENT OF 'K» WITHOUT lONIZ. 

11 REAL PART OF N-SQUARE (APRROX.) 

12 IMAGINARY PART OF N-SQUARED (APPROX.) 

13 REAL PART OF PROPAGATION CONSTANT (APPROX.) 

14 IMAGINARY PART OF PROP. CONSTANT (APPROX.) 

15 ATTENUATION CONSTANT WITH IONIZATION 

16 ATTENUATION CONSTANT WITHOUT IONIZATION 

17 CHANGE IN ATTENUATION DUE TO IONIZATION 

19 ■ ACOUSTIC SOURCE STRENGTH 



SUBROUTINE INGARO(X) 

IMPLICIT REAL*8 (D) 

IMPLICIT COMPLEXES (0) 

RFAI.^S'A X( 5) . n ,12.1 OSS,D,DI 
INTEGER*4 KLIC/0000/ 

EQUIVALENCE (QY1,RQY1) 

COMMON /INf|ER/C( 201 ,£N( 100) , Yl( ICO) ,Y2( 100 ) ,ARG(1C) , 
’S'VAL(IO) ,KCC, J 
COMMON /8L0C1/0UT(20) 

NAMEL 1ST /LISTA/CN, CL, BULK, WALL, LOSS 

NAMELIST /LI STB/SI GMA , T AU , W I N , WN I , W EN , WNE , GN , 11,12 

NAMELIST /LISTC/GN, D,DI ,AA,B1,B2,B 

NAMELIST /L I STD /WE, W I ,CCP, CLI I , CL I , WE I , W i E , C E , T AU I 
NAMELIST /LISTF/QA,QB,QC,QD,0E,QF,QG,QH,QI 
NAMELIST /LISTF/QYl ,QY2,QY3,QV 
NAMELIST /LISTG/QVI ,QVE 
NAMELIST /LISTH/RM,RO, THETA 



CN=C(11*SQRT(X(4) ) 

CL=C(2)/X(2) 

BULK=1. 54919*CL/CN^X( 1 )*=i=2 
WALL=2,0744434SQRT( X(1 )*CL*CN)/C(1C) 
LOSS=BULK+WALL 
IF(KCC.EQ.l) WRITE(6,LISTA) 



ENTER TABLE OF ENTRIES WITH VALUE OF ELECTRON TEMP 
AND OBTAIN VALUE OF Y2 . USE TO CALCULATE TAU. 



CALL ATSM( X( 5) , E N , Y2 , 100 , 1 , ARG , VAL , 5) 
CALL AL I ( X( 5) , ARG, VAL, I2,5,.1E-2,IER) 
SIGMA=C(3)=i'I2 

TAU=C(4)*X( 5)=A^^1.5-^=X(3 )4'SIGMA/CN*^2 



ENTER TABLE OF ENTRIES WITH VALUE OF ELECTrON TE‘'"> 



l6l 



no oo no oo oo oo on no 



c 



WEN 



AND OBTAIN VALUE OF Yl. USE TO CALCULATE 
CALL ATSM(X(5) , EN , Y 1 , 100 , 1 , ARC , VAL i 5 ) 
call ALI (X(5), ARG,VAL, 1 1 , 5t . lE-2, IER) 
WEN = C( 6)'-^X( 2)=i=SQRT( X( 5) )^I 1 
WNE=X(3 )/X(2)=i'C(7)^WEN 
WIN=C ( 5 )>!^CN/CL 
WNI=WIN»^X(3)/X(2) 

GN=V/NE+WNI 

IF(KCC.EO.l) WRITE (6,LISTB) 



1F(J.EQ.2) GO TO 401 



WE=C( 11 )*SQRT( X(3) ) 
WI=WE*SQRT(C(7 ) ) 



500 



RM=C(12)/X(4) 

RD=C( 13 ) /WE*SQRT( X( 5) ) 
THETA=R0/RM 

IFITHETAcGE.l. ) GO TO 5C0 

IF(KLIC.LE.2) WRITE(6,L ISTH) 

THETA=1 o 1 

KLIC=KLIC+1 

CONTINUE 

IF(KCC.EO.l) WRITE(6,LISTH) 



CLI I=1*/(X(3)^RM«*2*AL0G(THETA) ) 

CLI=1./ (l./CL! I+l./CL) 

IF(CL loLE.G.O) WRITE(6,LISTD) 

TAUI^l. 5492=^^00 I/CN’«=X( 1 )^*2+2.07444*SQRT( XU )=«^CLI*CN) 
=*=/C( 10 ) 

CCP=C( 15)#R0**3YX( 3) 

WE 1 = 3 c 6 2 E- 1 P^'-'X ( 3 ) * AL OG ( CCP ) /X ( 5 ) '>« 1 . 5 

WIE=C(7 )*WEI 

CE=C( 14 )-S0RT( X( 5) ) 

IF(KCC.EO.l) WRITE(6,LISTD) 



C 0 =(n.,l.) 

CA = U+QQ=^( LOSS + GN) /X( 1 ) 

QB=-QQ^VWNI/X( 1 ) 

QC=(TAU-QQ^WNE ) /X( 1 ) 

QD=-0Q=^WIN/X( 1 ) 

QE = -(WI/X( 1 ) )>!'^^2 + Q0*(TAUI+WIN+WlE)/X( 1) 

QF=(WI/X(1) )=}^=.'=2-QQ-Wl E/X( 1) 

QG=-0Q'7WEN/X( 1 ) 

0H=(WE/X(1 ) )=’^=:‘2-QQ + WEI /X( 1) 

0I=-(WE/X(1 ) + (8)^'(CE/CN)*=f'2 + Q0-(WEN+WEI )/X( 1 ) 

IF(KCC.EO.l) WRITE(6,LISTE) 



QYl = QE’i'QI-QH*QF 

ROY 1 = - ( ( W I E/ CN ) ❖« 2-0 . 6+ ( T AU I +W I N ) WEN+ WE I ) -W I E -W EN ) 
’i'/X( 1 )'i'--!'2 

QY2 = QC'-l'( Q0«0H-QE*QG) 

QY3=QB’!'(QGY0F-gD-0I ) 

QY=0A+( QY2+QY3 ) /QYl 
IF(KCC.EQ.I) write (6,LISTF) 



GUT(1)=REAL(0Y) 

OUT (2 )=AIMAG(QY ) 
QS1=CSQRT( QY) 
0UT(3)=REAL(QS1 )*X( D/CN 
0UT(4)=AIMAG(QS1 )*X( 1 ) /CN 



C 



IF(J.EO.l) GO TO 999 
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AOl CONTINUE 

D = X( ] )’i^X(3)^'=X(5>/C(8)/X(2)/X(A)/GN 
DI = 1. /( lo+D’f^D) 

AA=(GN+TAU)*D=i‘OI/X( 1 ) 

B1=GN*0*D*DI 

B2=TAU*DI 

B=( L0SS+B1-B2) /X(l ) 

IF(KCC.EQoi) WRITE(6,LISTC) 

OUT( 11 > = 1.-AA 
OUT(] 2)=B 

QX=CMPLX(0UT{1 1 ) ,0UT(12)) 
QS1=CSQRT( QX) 
0UT(13)=REAL(QS1)*X(1)/CN 
OUT( 14) =AIMAG( QSl ) *X( 1 ) /CN 



IF(J„NE.3) GO TO 999 



QVI=“( QD*OI-QG^-QF ) /QYl 
0UT{5)=CABS(QV1 ) 

T1=REAL( CVI ) 

T2=AIMAG(QVI) 

0UT(6) = ATAN2(T2,T1 ) 

QVE=~( OE-'‘QG-QD*QH) /QYl 
0UT(7)=CABS(QVF ) 
IF(KCC.EQ,1) WRITE(6,LISTG) 
T1=REAL (OVE) 

T2=AIY,AG(QVE ) 
0UT{8)=ATAN2(T2,T1 ) 



Fl=l e0 3 2 8^CL/CN*X( 1 )**2+l. 24467=I'SQRT ( X ( 1 )*CL*CN) /C( 10) 
F2=l .291'^CL/CN=i=X( 1 ) =^:*2 + ?„0744=f^S0RT( X ( 1 ) ’'^CL^i'CN ) /C ( lU ) 
QN1 = 1.+Q0>;'F1/X( 1 ) 

CN2 = l.-o4*0Q*F2/X( 1 )/ ( 1» + QQ«F2/X( 1 ) ) 

QN=QN1/QM2 

QS1=CSQRT(QN) 

0UT( 10)=AIMAG( QSl )*X(1 )/CN 



D0UT=0UT(4) /OUT ( 10) 

D0UT2=0G'JT-loD0 

0UT(9)=DCiUT2 

0UT( 15)=0UT(4) /SQRT(X{ 1 )/6.283183) 
0UTn6)=0UT( 10 )/SQRT( X( l)/6. 28 3183) 
0UT( 17)=0UT(4)-0UT( 1C) 

0UT( 19 )=CN>>-*2=i^TAU 



999 RETURN 
END 
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